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PROGRAM  SUMMARY 


This  report  presents  the  results  of  a  research  program  designed  to  investigate  the  generation 
of  ultraviolet  radiation  arising  from  the  non-equilibrium  thermodynamic  conditions  behind  strong 
shock  waves  in  air.  Such  conditions  can  arise  behind  the  bow  shock  of  vehicles  in  rapid  ascent 
through  the  atmosphere.  The  overall  objective  of  this  experimental  and  analytical  program  was  to 
support  the  development  of  a  capability  to  predict  reliably  the  non-equilibrium  ultraviolet  radiation 
from  such  flowfields. 

Specifically,  the  program  objectives  included: 

•  an  early  assessment  of  the  radiation  to  identify  the  radiating  species  and  to 
establish  its  dependence  on  velocity  and  altitude.  These  data  were  to 
support  instrument  design  issues  for  a  proposed  sounding  rocket  flight  to 
study  these  boost-phase  signatures. 

•  the  development  of  a  comprehensive  base  of  quantitative  spectral  radiance 
data  to  support  the  development  of  predictive  capabilities  by  the  boost-phase 
signature  community. 

•  the  development  of  a  model  to  describe  the  electronic  state  excitation  kinetics 
that  govern  the  nonequilibrium  ultraviolet  radiation  from  shock-heated  air. 

The  experiments  made  use  of  a  shock  tube  facility  and  spectro-radiometric  instrumentation 
available  at  Calspan.  The  system  was  capable  of  producing  shocked  air  conditions  over  the  actual 
range  of  velocity  and  altitude  of  interest  to  the  flight  program.  Both  spectrometric  and  radiometric 
data  were  obtained  as  the  shock  waves  swept  past  the  viewing  station  in  the  shock  tube.  Species 
radiating  in  the  ultraviolet  as  well  as  in  the  infrared  were  investigated  in  terms  of  spectral  intensity 
and  temporal  behavior.  Additionally,  shocks  in  various  nitrogen/oxygen  mixtures  were  used  to 
provide  a  broad  matrix  of  data  to  support  the  chemical  and  electronic  state  kinetics  developments. 

Following  is  a  listing  of  the  principal  results  that  have  accrued. 

1 .  Highly  resolved  spectra  of  the  nonequilibrium  UV  radiation  between  200  -  400  nm 
have  shown  the  NO  Gamma-band  radiation  to  dominate  over  a  smaller  contribution  from  the  O2 
Schumann-Runge  system. 


2.  The  dependence  of  this  nonequilibrium  UV  radiation  on  flight  vehicle  velocity  and 
altitude  was  established  by  direct  measurements  in  the  laboratory.  These  data  were  successfully 
applied  to  subsequent  flight  instrument  settings. 

3 .  Spectra  obtained  from  tests  using  pure  O2  permitted  the  contribution  of  O2  to  the  air 
spectrum  to  be  quantitatively  determined. 

4 .  Similarly,  the  temporal  profiles  in  pure  O2  recorded  by  the  radiometers  were  used  to 
obtain,  by  subtraction,  the  kinetic  behavior  of  the  dominant  NO-Gamma  band  radiation  behind  the 
shock  waves. 

5.  An  extensive  data  base  of  UV  nonequilibrium  spectral  radiance  profiles  has  been 
generated  for  use  by  the  signature  modelling  community  in  this  area.  Parameters  in  this  data  base 
include  shock  velocity,  pre-shock  pressure  (altitude)  and  test  gas  composition;  air,  pure  O2  and 
various  O2-N2  mixtures. 

6.  A  supporting  infrared  data  set  has  been  obtained,  which  provides  temporal  profiles 
of  the  global  NO  formation  kinetics  behind  shock  waves  in  air  and  O2-N2  mixtures. 

7 .  It  is  concluded  that  the  currently  accepted  rates  for  air  ground-state  chemistry  do  not 
predict  the  NO  generation  correctly. 

8.  A  simplified  excitation  model  is  suggested  to  describe  the  ultraviolet  radiation 
behind  the  shock  wave. 
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1.0  INTRODUCTION 


The  motivation  for  this  program  centered  on  the  development  of  capability  for  predicting 
the  ultraviolet  radiation  from  the  bow  shock  of  a  vehicle  ascending  between  40-60  km  at  velocities 
between  3-4  km/s.  The  objectives  of  the  program  included  providing  quantitative  measurements  in 
support  of  a  sounding  rocket  flight  experiment  the  identification  of  the  radiating  species,  and  the 
dependence  of  the  radiation  on  altitude  and  velocity.  In  addition,  basic  spectral  radiance  data  were 
acquired  in  support  of  modelling  this  nonequilibrium  radiation. 

The  non-equilibrium  aspect  of  the  radiative  signature  arises  from  the  fact  that  at  higher 
altitudes  the  decreased  density  and  collision  frequency  of  molecules  behind  the  shock  front  slows 
•both  the  chemical  reactions  that  establish  species  concentrations  and  those  that  govern  equipartition 
among  internal  energy  states:  electronic,  vibrational  and  rotational.  Thus,  under  some  conditions 
there  exist  strong  gradients  behind  the  shock  front,  consisting  of  species  concentrations  and  excited 
state  populations  far  from  equilibrium  with  the  local  (and  changing)  translational  temperature.  To 
describe  this  nonequilibrium  and  to  predict  its  dependence  on  velocity  and  altitude  requires  that  the 
relevant  kinetic  rates  and  mechanisms  be  known.  Thus,  a  major  goal  was  to  provide  a  model  to 
describe  the  electronic  state  kinetics  that  govern  the  nonequilibrium  ultraviolet  radiation  from 
shock-heated  air. 

The  experimental  aspects  of  the  program  centered  upon  the  use  of  a  shock  tube  and 
associated  spectrometric  and  radiometric  instrumentation,  whereby  radiative  signatures  from  shock 
heated  air,  O2  and  N2-O2  mixtures  were  recorded.  The  flight  conditions  of  interest,  namely,  shock 
velocities  of  3-4  km/s  at  altitudes  above  40  km  (2.25  torr)  can  be  exactly  duplicated  in  a  shock 
tube.  Thus,  the  measurement  of  the  radiation  profile  as  the  shock-heated  air  passes  a  given  station 
can  be  readily  and  simultaneously  achieved  for  a  number  of  wavelengths  and  with  different 
instruments.  This  procedure  was  the  basis  for  the  measurements  conducted  under  this  program. 

The  results  of  the  first  two  years'  work  are  embodied  in  two  AIAA  papers.  Refs.  1  and  2, 
which,  for  completeness,  are  included  in  their  entirety  as  Appendices  A  and  B.  The  shock  tube 
facility,  instrumentation,  calibration,  typical  data  and  their  description  and  analyses  are  all 
presented  in  these  papers. 

The  body  of  this  report  addresses  principally  the  efforts  of  the  subsequent  work.  These 
included  the  .search  for  additional  diagnostic  ob.servables  and  the  resultant  infrared  measurements 
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of  the  5.3  urn  NO  fundamental  vibronic  radiance  profiles.  Finally,  an  extensive  discussion  of  the 
gasdynamic  and  electronic  state  kinetic  modelling  and  analyses  is  presented. 


2.0  EXPERIMENTAL  CONSIDERATIONS 

2. 1  Search  for  Additional  Observables 


As  will  be  addressed  further  in  later  sections,  a  major  effort  was  undertaken  to  formulate  a 
model  to  describe  the  kinetics  of  the  population  of  the  electronically  excited  NO(A)  state.  This  state 
gives  rise  to  the  Gamma  band  system  which  dominates  the  nonequilibrium  UV  radiation  from 
shock-heated  air.  Because  the  model  contains  a  number  of  coupled  reactions  and  interacting 
species,  a  search  W2is  undertaken  to  determine  whether  any  other  relevant  species  gave  rise  to 
observables  that  could  be  used  for  additional  kinetic  diagnostics.  Two  were  identified:  atomic 
oxygen  and  the  electronically  excited  N2(A)  state.  Each  of  these  playja  key  role  in  the 
nonequilibrium  radiation  overshoot  kinetics;  serious  effort  was  expended  to  evaluate  their  potential 
as  experimental  observables. 

The  technique  of  atomic  oxygen  absorption  near  130  nm  in  the  vacuum  ultraviolet  is 
actively  being  used  in  shock  tube  kinetic  studies,  as  reported  in  the  1989  Shock  Tube 
Symposium.  Details  of  this  technique,  especially  those  of  the  absorption  lamp  source  were 
obtained  and  found  to  be  tractable.  However,  final  calculations  showed  that  the  concentration  of 
O-atoms  in  the  present  experiments  was  much  too  high  for  these  absorption  measurements. 
Atomic  O  absorption  would  correspond  to  an  optically  thick  gas,  which  precluded  its  use  ais  a 
kinetic  diagnostic,  and  resulted  in  abandoning  this  effort. 

The  use  of  laser-induced  fluorescence  techniques  to  obtain  instantaneous  concentrations  of 
the  excited  N2(A)  state  behind  shock  waves  was  also  considered.  The  experimental  plan  would  be 
to  absorb  the  laser  beam  by  transitions  from  the  A  to  B  state  and  to  observe  subsequent  fluorescent 
emission  to  other  vibrational  A  states.  These  bands  are  in  the  N2  First-Positive  B%  -  A^£  system, 
between  600  -  800  nm. 

This  design  effort  involved  a  detailed  study  of  the  available  energy  levels  of  these  states, 
the  transitions  in  various  wavelength  regions  and  the  resultant  line  der  ities  per  wavelength 
interval.  Overlap  integrals  were  taken  into  account  to  optimize  fluorescent  signals  while 
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minimizing  laser  power  requirements.  The  use  of  cxcimer  and  Y AG-pumped  dye  lasers  was 
evaluated,  as  either  could  be  procured  for  the  term  of  the  experiments. 

The  results  of  these  analyses  showed  that  this  measurement  was  both  high  risk  and  costly, 
and  that  its  implementation  could  adversely  affect  the  balance  of  this  program.  Several  factors 
ruled  against  this  measurement.  A  major  one  was  that  the  temperature  of  the  gas  during  the  region 
of  interest  was  high  (>  4000  “K),  resulting  in  a  significant  thermal  Boltzmann  population  of  both  B- 
and  A-state  levels.  Further,  because  they  are  separated  by  only  1.2  ev  (B  =  7.4  ev,  A  =  6.2  ev), 
their  populations  are  not  much  different,  typically  10^  vs  10^  per  cc.  Thus,  when  attempting  to 
measure  the  fluorescent  component  over  a  manifold  of  transitions,  the  inability  to  differentiate 
against  spontaneous,  thermal  transitions  becomes  prohibitive.  Efforts  on  this  task  were 
terminated,  and  further  diagnostic  efforts  were  directed  toward  the  improved  infrared 
measurements  described  below. 


2.2  Infrared  Measurements  of  NO 


From  concurrent  analyses  it  became  apparent  that  the  kinetics  of  the  formation  of  NO 
behind  the  shock  played  a  critical  role  in  the  overall  scheme  for  the  excitation  of  the  electronic 
NO(A)  state.  Earlier  in  the  program, 2  preliminary  measurements  of  the  vibronic  bands  of  NO  at 
5.3  pm  had  shown  that  this  system  could  be  used  to  monitor  global  (ground  state)  NO  kinetics. 
These  early  profiles  were  measured  with  a  standard  IR  radiometer  comprising  a  HgCdTe  detector, 
bandpeiss  filter  and  slits  for  FOV  definition.  Figure  1 1 ,  page  B-9  depicts  several  such  profiles,  and 
illustrates  the  complication  of  a  ’’foot”  of  signal  preceding  the  arrival  of  the  shock  wave  at  the 
observation  port  FOV. 

This  foot  generally  arises  from  scattered  light  from  the  highly  luminous  shock  front,  as  it 
approaches  the  measurement  station.  Scattering  can  occur  from  any  defect  (scratch,  chip,  film)  on 
either  of  the  shock  tube  windows,  or  from  any  slits,  baffles  or  radiometer  tube  walls  that  are 
illuminated  by  the  shocked  gas  front  and  which  are  ’’seen”  by  the  detector.  Even  second  reflections 
can  give  rise  to  the  scattered  component.  It  should  also  be  noted  that  after  the  shock  peisses 
through  the  FOV,  any  candidate  scattering  site  or  surface  is  doubly  illuminated,  giving  rise  to  a 
larger  scattered  signal  component,  but  which  is  now  buried  in  the  true  signal.  The  correction 
would  be  at  least  equal  to  twice  that  of  the  foot  seen  when  the  shocked  gas  is  not  in  the  FOV. 
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A  new  radiometer  was  designed  to  surmount  these  difficulties.  A  schematic  diagram  of  the 
optics  is  shown  in  Figure  1.  A  major  feature  is  that  the  entire  radiometer  is  separated  from  the 
shock  tube.  The  physical  apenures  that  define  the  very  precise  FOV  are  each  imaged  onto  the 
shock  tube  windows  by  a  CaF2  lens.  One  stop  is  the  detector  chip  itself,  and  the  second  is  a  slit 
made  of  razor  blades,  well  removed  from  illumination  by  the  advancing  shock  wave.  It  is  sized 
and  positioned  so  that  its  image  on  the  lower  window  is  appropriate  for  the  desired  FOV.  In  this 
way  there  are  no  physical  sites  from  which  scattering  can  occur.  The  windows  were  removed  after 
each  test  to  preclude  scratching  while  the  tube  was  being  scrubbed.  Implementing  these  features 
resulted  in  a  dramatic  reduction  in  scattered  radiation,  as  evidenced  by  the  lack  of  the  pre-shock 
"foot"  in  the  temporal  radiance  profiles.  Figure  2  illustrates  the  improvement. 

An  InSb  detector  was  used  in  the  new  radiometer.  It  was  found  that  even  though  the  long 
wave  roll-off  of  the  detector  response  clips  one  wing  of  the  NO  band,  the  higher  detectivity  results 
in  a  net  gain  in  SNR.  The  combined  filter-detector  bandpass  function,  between  5  -  5.5  |im  was 
determined  by  separate  bench  experiments.  A  standard  blackbody  source  through  a  scanning 
monochromator  was  calibrated  and  then  used  to  establish  the  wavelength  dependence  of  the 
radiometer,  as  shown  in  Figure  3. 

System  calibration  was  achieved  by  means  of  a  chopped,  standard  blackbody  which 
completely  filled  the  FOV  and  whose  intensity  spanned  the  range  of  those  for  the  test  gases. 
Figure  4  presents  a  plot  of  the  blackbody  radiance  convolved  with  the  filter  function  and  the 
detector  response  function  (jlBBFxD>,d>.,  wcm'^sr)  against  the  recorded  signal  levels  in  millivolts. 

System  linearity  is  evident,  and  these  data  were  used  to  obtain  an  overall  calibration  constant  that 
was  applied  to  all  subsequent  profile  measurements.  By  these  means  the  in-band  radiance  from  the 
NO  (formed  and  heated  behind  the  shock  wave)  can  be  directly  related  to  that  from  a  blackbody  of 
known  temperature.  The  analysis  of  these  data  and  their  application  to  the  kinetic  modelling  are 
discussed  in  a  later  section. 


3.0  ANALYTICAL  AND  MODELLING  CXINSIDERATIONS 
3.1  Equilibrium  Radiation  of  NO  and  O2 

The  time-dependent  measurements  of  the  ultraviolet  radia*'on,  as  described  in  Appendix  A 
page  3.  can  be  used  to  obtain  an  equilibrium  radiation  value  at  long  times.  These  can  then  be 
compared  with  previous  equilibrium  experiments  and  calculations  by  reducing  the  results  to  an 
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evaluation  of  the  electronic  transition  moments  for  the  NO  Gamma  and  O2  Schumann-Runge 
systems. 


To  accomplish  this,  the  radiation  from  a  single  spectral  vibrational  band  of  an  electronic 
band  system  is  written  as 


^  ) 


(1) 


where 


3  <2.  2 

l(^  TT  ce  gc  10  =  3,1C>‘10 


4  / 

W  CAY}  / SfU 


6d  = 

S’lRel^ 


Q. 


I 


number  density  of  molecules  (cm‘-^) 
fraction  of  molecules  in  upper  electronic  state 
wavenumber  of  radiation  (cm"') 

sum  of  squared  dipole  moments  for  transitions  from  degenerate  upper 

states,  atomic  units 

degeneracy  of  upper  state 

vibrational  energy  of  upper  state 

vibrational  partition  function  of  upper  state 


By  representing  each  vibrational  transition  (v'  — »  v")  as  a  simple  Q  branch,  and  assuming  that 
IR^I^  can  be  taken  as  a  constant,  the  spectral  radiance  can  be  approximated  as 


1  Be 

l^e' 

<  1 

f 

\  J  W 

IkT 

CJnr^a^  CrrT' 

(2) 


O  f-on,  co^CJyi^" 


co>Cc) 


vV 


where 


Be',  6.  =  rotational  constants  for  upper  and  lower  states 

^vV"  ~  bandhead  wavenumber  for  v',  v'  transition 

If  Be’>  Be  ",  the  equation  applies  forco^c^^y,  and  ®  for  .  The  dependence  of 

B  on  v  is  nof  included  in  the  approximation.  The  total  spectral  radiance  is  then 
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(3) 


VV" 

The  spectral  radiance  calculated  with  this  approximation  is  shown  in  Appendix  A,  Figures  12  and 
13  for  the  NO  Gamma  and  O2  Schumann-Runge  systems. 

The  in-band  radiation  that  should  be  measured  by  the  calibrated  detector  for  a  single  band 
system  is  then 

S  =  X  1  J  Cui)  ^  (co)  D  Cu}')  dco  W  (4) 

where  £  is  the  diameter  of  the  shock  tube,  F(6j)  is  the  filter  transmission  function  and  D  (to)  is  the 
detector  response  function.  F(cj)  for  the  230  and  260  nm  filters  is  shown  in  Figure  5  of  Appendix 
B.  The  detector  that  was  used  for  the  broad  230  nm  filter  has  maximum  sensitivity  at  205  nm. 
The  sensitivity  decreases  by  20%  at  185  nm  and  at  250  nm.  At  longer  wavelengths  the  sensitivity 
drops  exponentially,  decreasing  by  a  factor  of  10  at  300  nm.  The  sensitivity  of  the  detector  that 
was  used  for  the  narrow-band  260  nm  mezisurements  was  taken  as  constant  over  the  band-pass. 

The  integrals  in  Equation  (4)  were  calculated  as  functions  of  temperature  for  the  NO 
Gamma  and  O2  Schumann-Runge  systems.  The  absolute  calibration  was  obtained  by  calculating 
similar  integrals  for  a  black-body  source  of  radiation  and  following  the  calibration  procedure 
described  in  Section  2.2. 

This  results  in  calculated  functions  of  temperature,  MnO  and  Mo2>  where  these  are  related 
to  SnO  and  So2  by 

Swo  ~  ^No 

From  Equations  (4)  and  (5)  it  follows  that  the  in-'  md  radiation  to  be  expected  from  the  gas 
mixture  would  be 
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Smos  =  =n,^l'^eUo  i  (R«.lei  I  IMoj. 


(6) 


or 


>MIX 


>IoiMo2.X 


=  iReloa  ^  l*^e,l#s/o  (  y^t^o  \ 

Hoi  Vto2  / 


(7) 


Since  MNoand  M02  are  known  from  the  calculation,  and  nNoand  no2  are  known  from  the 
equilibrium  shock  calculation,  the  experimental  value  of  Smix  divided  by  yvo2Mo2^  caa  be 
plotted  vs  (Mno  nN0^02n02)-  intercept  and  slope  of  the  straight  line  yield  the  |ReP  for  O2 

and  NO.  These  plots  are  shown  in  Figures  5  and  6  for  the  230  nm  and  260  nm  results  for  O2  and 
NO.  Measured  vs  calculated  values  of  the  equilibrium  in-band  radiance  for  230  and  260  nm  are 
plotted  in  Figures  7  and  8  of  Appendix  B,  using  the  derived  values  of  IR^P.  (The  notation  of 
Appendix  B  is  slightly  different  from  that  used  here.) 


3.2  Separation  of  Noneauilibrium  0->  and  NO  radiation 

The  equilibrium  ultraviolet  radiation  from  the  shock  tube  is  a  well  understood  mixture  of 
O2  Schuniann-Runge  radiation  and  NO  Gamma  radiation,  as  discussed  in  Section  3.1.  The  5%  O2 
experiments  produce  practically  no  O2  radiation,  the  22%  O2  tests  provide  about  equal  amounts  of 
NO  and  O2  contributions,  and  the  40%  O2  tests  provide  appreciably  more  Schumann-Runge  than 
Gamma  contributions.  In  the  nonequilibrium  region,  the  NO  radiation  is  the  dominant  contributor 
in  almost  all  cases.  However,  it  is  necessary  to  make  some  correction  to  the  total  radiation 
measurements  for  the  O2  content  in  order  to  obtain  the  NO  overshoot  history. 

Experiments  with  100%  oxygen  show  that  the  O2  radiation  does  not  overshoot  its 
equilibrium  value,  but  rises  smoothly  to  a  final  plateau  (see  Figure  5  Appendix  A).  This  feature 
was  observed  in  work  at  AVCO^  in  the  late  1950’s,  and  attributed  to  the  rapid  predissociation  of  the 
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oxygen  excited  B  state.  Thus  the  B  state  remains  in  local  thermodynamic  equilibrium  with  the 
oxygen  atom  population,  so  that  the  number  of  O2  molecules'cm^  in  the  B  state  is  given  by 


where  K  is  the  equilibrium  constant  for  the  reaction  and  the  last  factor  in  the  equation  is  an 
approximate  representation  (appropriate  over  the  range  3000  to  7000  K)  of  the  required  correction 
to  the  rigid-rotator  simple-harmonic-oscillator  model.  The  O2  radiation  is  taken  to  be  proportional 
to  no2(B)7  as  determined  from  the  equilibrium  measurements. 

We  have  assumed  that  the  O2  S-R  system  follows  the  same  behavior  in  the  radiation  from 
N2/02  mixtures,  with  the  02(B)  state  staying  in  local  equilibrium  with  the  oxygen  atom 
concentration  at  the  translational  temperature.  The  oxygen  atom  concentration  behind  the  shock 
was  calculated  from  a  normal-shock  nonequilibrium  code^,  and  this  population  was  used  to 
determine  the  B  state  participation  in  the  observed  radiation.  This  amount  of  radiation  was 
subtracted  from  the  total  observed,  and  the  remainder  was  attributed  to  the  NO  Gamma  system. 

The  results  of  this  subtraction  are  shown  in  Figure  6  of  Appendix  B  for  a  matrix  of  nine 
runs.  The  total  radiation  and  the  portion  attributed  to  NO  Gamma  are  shown.  In  the  5%  O2 
experiments  the  corrections  are  negligible,  and  for  the  22.3%  O2  tests  the  corrections  are  very 
small  in  the  nonequilibrium  region,  except  at  the  slowest  shock  velocities.  For  the  40%  cases  at 
low  velocities  there  is  an  appreciable  correction. 


3.3  Analysis  of  Nonequilibrium  Infrared  Measurements 

The  nitric  oxide  5.3  um  radiation  intensity  is  expected  to  be  proportional  to  the  NO 
concentration  and  a  relatively  weak  function  of  the  temperature*.  Thus  the  measurement  of  this 
radiation  can  be  used  to  determine  the  time-history  of  the  NO  concentration  behind  the  shock  wave. 
This  time-history  of  the  ground-state  chemistry  behind  the  shock  wave  is  a  critical  input  to  any 
calculation  of  the  anticipated  time  history  of  the  ultra-violet  radiation.  The  infra-red  measurements 
obtained  during  the  first  phase  of  this  prog'tm  demonstrated  the  ability  to  determine  NO 

*  The  NO  is  optically  thin,  and  it  is  show-n  in  Reference  3  that  NO  vibrational  relaxation  does  not  appreciably  affect 
the  interpretation  of  the  measured  NO-IR  radiation.  The  temperature  dependence  is  given  in  Appendix  D. 
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concentration,  but  did  not  contain  enough  resolution  to  resolve  this  time  histor>'.  In  the  second 
phase  of  the  program  the  sensitivity  of  the  infrared  detection  system  was  improved  to  give  a  better 
measure  of  this  ground-state  NO  chemistry,  to  use  as  a  guide  in  interpretation  of  the  ultraviolet 
data. 


A  very  general  picture  of  the  observed  infra-red  radiation  is  given  in  Figure  7  showing 
plots  of  observed  intensity  vs  lab  time  behind  the  shock.  (Full-page  copies  of  these  graphs  are 
given  in  Appendix  C).  The  three  columns  correspond  to  the  three  O2/N2  mixture  ratios  that  were 
investigated,  and  the  results  are  arranged  vertically  according  to  velocity.  The  time  histories  show 
a  rise  to  maximum,  an  approach  to  equilibrium,  and  the  disturbance  caused  by  arrival  of  the 
hydrogen/oxygen  interface,  resulting  in  OH  and  H2O  radiation. 

The  radiation  level  attained  at  equilibrium  is  proponional  to  the  NO  concentration  and  a 
weak  function  of  the  equilibrium  temperature.  Both  the  NO  concentration  and  the  temperature  can 
be  calculated  if  the  ground-state  chemical  rates  are  known.  It  is  shown  in  Appendix  D  that  the 
temperature  dependence  that  is  expected  for  NO  radiation  measured  with  the  filter-detector  system 
that  was  used  is: 


f  (T)  =3.SlXlo'^j^l  +•3.2.1x10^(1—3000^ 


(9) 


This  temperature  dependence  can  be  confirmed  by  dividing  the  measured  equilibrium  intensities  by 
f  and  plotting  the  results  vs  the  NO  concentrations.  As  shown  in  Figure  1  of  Appendix  D,  this 
demonstrates  the  proper  linear  dependence  on  NO  concentration,  and  provides  a  relation  that  can  be 
used  to  interpret  the  nonequilibrium  data. 


Also  shown  in  Figure  7  are  calculations  of  the  expected  .NO  IR  radiation.  These  calculated 
results  were  obtained  by  utilizing  a  7-species.  lO-reaction  model  in  a  normal-shock  kinetics  code 
to  determine  the  time  history  of  the  NO  concentration  and  temperature  behind  the  shock.  These 
values  were  then  used  with  Equation  (12)  of  Appendix  D  to  determine  the  radiation.  The  rates  that 
were  used  in  these  calculations  are  those  suggested  by  Wray^,  as  corrected  by  Camac  et  aF.  The 
Camac  corrections  are  the  deletion  of  the  reaction  N2+02^2N0.  and  a  decrease  of  a  factor  of  2.5 
in  the  rate  for  .\2+0—N0+N.  The  Camac  rates  are  expressed  as 


k. 


(10) 
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and  are  given  in  Table  1.  These  rates  were  considered  to  be  the  most  relevant  ones  available  for  the 
density  and  shock-velocity  range  of  interest  in  the  present  experiments.* 


jjj.  =  AT^g-C/T  cm^/molecule-sec  or  cm^/molecule-sec 


A 

B 

C 

(°K) 

1. 

O2  -t-  0  ^  20  0 

1.5  X  10-4 

-1.0 

59,380 

2. 

O2  -<■  O2  20  O2 

5.4  X  10-5 

-1.0 

(( 

3. 

ISHSBHRim 

1.2  X  10-5 

-1.0 

n 

IH 

6.8  X  10-2 

-1.5 

113,200 

5. 

N2  +  N2  2N  +  N2 

8.0  X  10-7 

-0.5 

/I 

6. 

N2  +  M  ^  2N  -(-  M 

3.2  X  10-7 

-0.5 

ft 

n 

NO-t-M^N-t-O-t-M 

1.7  X  10-3 

-1.5 

75,490 

8. 

N  -t-  O2  -  NO  -1-  0 

2.2  X  10-14 

H-l.O 

3,565 

9. 

5  X  10-11 

0 

38,016 

10. 

NO-*-  -t-  e-  -  N  -t-  0 

9.0  X  10-3 

-1.5 

0 

TABLE  1.  CAMAC  RATES  USED  FOR  GROUND-STATE  CALCULATIONS 


A  recent  publication*  gives  an  ex’- erimental  result  for  k3  (with  M=N2)  of  3.4xl0l*  (t  25%)T-' 
exp(-59,380/T),  measured  from  24(X)  to  4100K.  This  is  atwut  half  the  Camac  value.  The  most 
recent  experiments  with  NO  formation^  suggest  that  k9  =  1.8  x  lO*"^  exp  (-38,370)/T,  about  6 
times  larger  than  Camac’s  value. 
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The  general  conclusion  that  can  be  drawn  from  comparison  of  calculations  and  experiments 
in  Figure  7  is  that  although  the  calculated  results  show  profiles  similar  to  the  experiments,  they  do 
not  predict  the  magnitude  or  position  of  the  overshoot  correctly.  The  calculations  fail  most 
noticeably  for  low  velocities  and  for  the  40%/60%  02i'N2  mixture.  Camac’s  results^  showed 
similar  problems,  but  those  could  not  be  resolved  because  of  possible  boundary-layer  effects  in  his 
shock  tube,  which  had  a  1.5-inch  diameter.  In  the  present  experiments,  with  a  3-inch  diameter 
tube,  boundary-layer  effects  do  not  contribute  to  determining  the  time  scale  behind  the  shock  at  the 
times  involved.  In  addition,  many  of  the  calculations  do  not  demonstrate  a  radiation  maximum, 
whereas  an  experimental  overshoot  is  clearly  shown. 

Because  of  the  importance  of  the  ground-state  chemistry  in  explaining  the  measured 
ultraviolet  radiation,  some  effort  was  made  to  adjust  rates  to  provide  better  agreement  with  the 
infra-red  experiments. 

As  pointed  out  by  Camac  et  aF,  the  general  features  of  the  time  history  of  the  IR  radiation 
can  be  characterized  by  the  three  temporal  parameters  shown  in  Figure  8.  These  are  an  incubation 
time,  a  time  to  half  maximum,  and  a  time  to  maximum.  The  experimental  values  for  these 
parameters  are  plotted  in  Figures  9, 10,  and  11.  The  ratio  of  peak-to-equilibrium  radiation  is  shown 
in  Figure  12.  Each  figure  shows  experimental  curves  for  the  three  mixtures,  as  well  as  the  values 
obtained  from  normal-shock  calculations  using  the  Camac  rates.  Calculated  profiles  for  the  40% 
€>2/60%  N2  mixture  did  not  show  maxima. 

INCUBATION  TIME 

The  early  formation  of  NO  behind  the  shock  is  governed  by  reaction  9,  which  depends  on 
oxygen  atoms  to  produce  the  NO.  (Reaction  9  also  produces  a  nitrogen  atom,  which  is 
immediately  converted  to  another  NO  molecule  leaving  a  replacement  for  the  O  atom,  all  through 
reaction  8.)  Thus  there  is  a  delay  in  the  NO  production  until  O  atoms  are  formed  by  reactions  2 
and  3.  These  dissociation  reactions  are  delayed*®  by  vibration-dissociation  coupling,  an  effect  that 
is  not  included  in  the  present  calculations.  Thus  a  complete  description  of  the  incubation  time  of 
NO  formation  requires  an  accurate  description  of  the  early  process  of  oxygen-atom  dissociation. 
The  need  for  a  longer  incubation  time  led  Camac  to  discount  the  N2+02^2NO  reaction,  which 
would  supply  NO  at  very  early  times  (he  also  had  other  experimental  reasons  for  removing  this 
reaction).  However,  it  is  clear  that  the  Camac  rates  do  not  provide  long  enough  incubation  times  at 
low  velocities. 
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TIME  TO  HALF-MAXIMUM 


Apart  from  the  relatively  short  incubation  time,  the  time  to  half-maximum  is  governed  by 
the  rate  of  production  of  NO  by  reaction  9.  It  was  on  the  basis  of  this  measurement  that  Camac  et 
al  suggested  a  reduction  of  a  factor  of  2.5  in  this  rate*.  In  their  examples,  they  reduced  the 
coefficient  from  Wray's  value  of  7x10*^  to  2.8x10*^,  as  indicated  in  Table  1.  This  correction 
provides  reasonable  agreement,  except  for  the  low-speed  shock  waves;  for  these  cases  the  Camac 
rates  predict  a  longer  time  to  half  maximum  than  shown  in  the  experiments,  and  do  not  predict  an 
overshoot. 

TIME-TO-PEAK 

Figures  11  and  12  describe  the  time-to-peak  and  the  ratio  of  peak-to-equilibrium  radiation. 
The  Camac  rates  do  not  provide  results  consistent  with  the  experiments  for  either  parameter.  In 
many  cases  the  calculated  intensities  do  not  show  a  maximum,  as  is  evident  in  Figure  7. 

By  changing  the  reaction  rates  we  have  obtained  calculations  in  reasonable  agreement  with 
the  experimental  results,  but  only  by  changes  that  are  inconsistent  with  known  values  of  these 
rates.  For  example,  by  changing  rates  of  ki,  k2,  k3,  and  k9.  the  results  shown  in  Figure  13  can  be 
obtained.  These  calculations  still  do  not  provide  an  accurate  fit  to  all  of  the  data,  and  require  use  of 
rates  that  are  in  complete  disagreement  with  accepted  rates.  However,  we  have  used  these 
calculations  in  Section  3.5  to  make  comparisons  with  ultraviolet  radiation  because  of  their  close 
match  to  the  infrared  results. 

The  possible  changes  in  rates  that  would  provide  results  consistent  with  all  of  the  data  have 
not  been  pursued  further  in  the  present  work.  Instead,  the  infrared  data  have  been  used  as  a 
guideline  in  interpreting  agreement  of  theory  and  experiment  in  the  ultraviolet  measurements. 

3.4  Excitation  Equations  for  the  A^  1+  State  of  NO 

The  mechanism  by  which  the  NO(A)  state  is  excited  behind  the  shock  wave  is  not  knowjt,^ 
but  it  is  certain  to  involve  a  complex  set  of  competing  reactions.  In  particular,  energy  transfer  from 


Camac  et  al  arrived  at  the  conclusion,  based  on  a  first-order  analysis,  that  the  same  effect  could  also  be 

accomplished  by  decreasing  the  rate  of  reactions  2&3  by  a  factor  of  2.5.  which  would  decrease  the  supply  of 
oxygen  atoms.  This,  they  arg_.d.  would  decrease  the  NO  production  rate.  Although  this  is  uue  at  the  very 
earliest  times,  we  have  found  that  decreasing  the  02  dissociation  rate  increases  the  temperature  at  any  given 
position  behind  the  shock,  and  the  net  effect  is  an  increase  in  the  NO  formation  rate,  which  is  very  temperature 
sensitive. 


the  metastable  N2(A)  state  is  usually  considered  as  a  parallel  path  for  NO(A)  excitation.  A  listing 
of  the  energies  involved  in  relevant  chemical  states  is  shown  in  Table  2. 


N  +  N 

1 13,243  “K 

N  +  O 

86,321 

NO(A)  +  NO(X) 

74,375 

N2(A) 

72,234 

NO(A) 

63,259 

0  +  0 

59,400 

N 

56,621 

0 

29,700 

NO  +  NO 

22,232 

NO 

11,116 

N2 

0 

02 

0 

TABLE  2.  ENERGIES  ASSOCIATED  WITH  CHEMICAL  SPECIES 


Mechanisms  for  the  high-temperature  excitation  of  the  NO(A)  state  have  been  examined  in 
References  11  and  12.  The  simplest  mechanisms  suggested  are  direct  excitation  by  electron 
collision  or  atom  collision,  recombination  of  N  and  O  to  form  NO(A),  and  recombination  of  two  N 
atoms  to  form  N2(A).  The  near-resonant  exchange  between  NO(A)  and  N2(A)  is  usually  assumed 
to  occur  very  rapidly  compared  with  other  reactions. 

A  representative  set  of  reactions  is; 


1) 

NO+e 

NO(A)+e 

2) 

N+O+M 

NO(A)+M 

3) 

N2+e 

N2(A)+e 

4) 

N2+N 

N2(A)+N 
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5)  N+N+M  N2(A)+M 

6)  N0(X)+N2(A)  N0(A)+N2(X) 

7)  NO(A)  ^  NO(X)+hV 

Reaction  7,  the  spontaneous  emission  of  radiation  from  the  NO(A)  state  is  the  signal  that  is 
observed.  The  radiative  lifetime  of  this  state  is  known  to  be  0.2  ps,  so  that  k7f  =  SxIO^sec*.  The 
reaction  is  written  in  a  single  direction  because  the  radiating  gas  is  optically  thin.  Reactions  2  and  S 
are  written  as  three-body  reactions,  but  two-body  reactions  could  also  contribute.  In  the 
nonequilibrium  region  these  reactions  proceed  in  the  reverse  direction,  as  quenching  reactions.  It 
should  be  emphasized  that  a  great  many  other  paths  of  excitation  are  possible,  many  including 
intermediate  states. 

The  rate  of  increase  of  the  number  density  of  NO(A)  and  N2(A)  molecules  can  be  written 
for  this  set  of  reactions  as 

dt 

Vl  NoYlt^xA  "  ^  No  A  (11) 

(12) 


A  great  simplification  in  the  investigation  of  these  equations  comes  from  the  fact  that 
Reactions  1  through  7  do  not  appreciably  affect  the  number  densities  of  the  unexcited  species  nN2, 
no2-  "NO.  tiN.  no,  n^ot  and  n^.  This  follows  from  the  fact  that  the  excited  species  number 
densities  nNOA  and  nN2A  are  about  e- ”^0000/7  5x10'^  to  5x10-10  of  thg  corresponding  ground 

state  species  number  densities.  To  produce  ground  state  species  at  an  appreciable  rate  from  these 
excited  states  would  require  rate  constants  that  were  proportionately  greater  than  the  rates  used  for 
ground-state  chemi'  ry.  Likewise,  the  chemistry  represented  by  these  reactions  does  not  affect  the 
enthalpy  or  temperature  history  of  the  gas  behind  the  shock  wave. 
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With  this  simplification  Equations  (11)  and  (12)  can  be  rewritten  as 

dt 

dfc  HZA--  [^3f 

r  /  T 

-  '^NZA  6  ID  "^MOA 

The  rate  constants  are  specified  as  functions  of  the  temperature,  and  the  species  number  densities 
inside  the  brackets  are  known  functions  of  the  distance  behind  the  shock  wave.  It  follows  that  if 
the  ground-state  chemistry  behind  the  shock  is  known,  then  for  a  given  choice  of  rate  constants, 
the  bracketed  terms  can  be  considered  as  known  functions  of  the  time,  and  the  two  equations  can 
be  integrated  numerically  to  determine  njvoA  and  nN2A- 


The  relation  between  the  forward  and  reverse  rates  can  be  written  explicitly  a: 

kio.  _ /yt<v/c  \  » /ativox  ^ <g  .  ksit  \  =  ^ 

kl-F  \dM0A/  ^  \clM2Ay 

/<iNOK\  p 

>tN_^  _/  Kk  \  ^  /  K>isa  \  _  / Tt, w  \  / c6twa.x\  p 

k.5-^  i^NiA/gj  \^'^24q 

where  the  terms  dj  are  the  degeneracy  of  the  i*  state: 


kniz 

k+f 


E~mza 

kT 


dN2X=l;  dN2A=3;  dNOX=^’  dNOA=2 

Since  very  little  is  known  about  the  values  of  the  rate  constants,  there  is  some  merit  in 
simplifying  the  equations  further.  It  is  first  observed  that  the  sum  of  Equations  (13)  and  (14) 
provides  an  equation  that  does  not  involve  the  exchange  reactions: 
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r 


iitr 


«t-  "  [^ifL^e'^^s.nP^t^  -f-  ^7fJ  ^KfoA 

"  [j^3a^e  ^  fJZA 

If  it  is  next  assumed  that  the  exchange  reactions  are  very  fast  compared  with  the  other  excitation 
reactions,  nNOA  and  nN2A  will  remain  in  a  local  thermal  equilibrium  with  each  other,  so  that 

£M2<k  ~E»jca  8ot6 

yiNOA  -  Q  -  ^  .  YImo  ^  =-i-'^kw£0 

V^LnzA  cLnzA  cLhiox  il'Wa-  ^  ^M2. 


Q  is  then  a  known  function  of  time  behind  the  shock  wave,  with  a  value  between  0  and  about  0.1. 


Substituting  Equations  (15)  and  (17)  in  Equation  (16)  results  in 


=  F-  C5.£ 
At 


(18) 


where  Z,  F,  and  G  are  defined  as 


^  ^NOA*^fslZA  > 


;  Yl  =- 

( <  +  a)  (  t  i-es> 


=  r  f-  ka.-P  ‘)lK/>le>lM 

L  Tz 

i-  ^,p  Km  Km  i-'icjl 

n  J 


(19) 


(20) 


(21) 
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and  pi ,  P2.  and  ps  are  defined  as 


h 


_ E^fs/OA 

/  fcT 

(Xjnoa.  ->1 

- O  .  Al/NO 

A^tJOX 


(22) 


_  ^NfOA 

,  *^ls/0  I  I  • 

cLi^ox  K/MO 


NO  A 


NOX 


e  .  ^No 


X/»/2.\  Vt^| 


(23) 


(24) 


Pi  is  the  value  of  nNOA  in  local  thermodynamic  equilibrium  with  nNOX>  P2  is  the  value  of  unoa  in 
LTE  with  nNno,  and  p3  is  the  value  of  unoa  in  LTE  with  niM^,  ail  under  the  condition  that  nNOA 
and  nN2A  are  in  local  thermodynamic  equilibrium  with  each  other.  In  final  equilibrium,  pi=p2~P3- 
Equation  (18)  provides  the  differential  equation  for  nNOA(0- 

A  final  simplification  that  is  appropriate  over  much  of  the  relaxation  zone  is  the  assumption 
that  Z  is  near  local  equilibrium  ••  i.e.,  dZ/dt  is  small  compared  with  the  terms  F  or  (Q/(l-t<^))GZ  in 
Equation  (18).  This  is  true  when  changes  in  the  ground-state  chemistry  are  slow  compared  with 
the  electronic  excitation  mechanisms.  In  this  case  Z  can  be  taken  as 

Z  =  Zle  =  (F/G)  ((1+Q)/Q),  or  (nNOA)LE  =  (F/G).  (25) 

This  local  equilibrium  value  is  somewhere  between  the  local  values  of  pi,  p2,  and  P3,  depending 
on  the  relative  values  of  the  numerators  of  the  three  terms  in  Equation  (21).  If  k?  has  an 
appreciable  value  compared  with  the  other  terms,  it  serves  to  lower  the  local  equilibrium  value  of 
tiNOA- 


19 


When  the  ground-state  chemistry  arrives  at  final  equilibrium,  pi=p2=P3  =  (nNOA)FE-  (see 
Equations  22-24),  and 

=  L  =  _f - -  (26) 

If  k7«F/(nNOA)FE-  then  (nNOA)LE=(nNOA)FE.and  the  equilibrium  radiation  is  not  "collision 
limited." 

3.5  Rate  Constants  for  Electronic  Excitation 

As  discussed  in  Section  3.4.  References  11  and  12  give  assessments  of  the  state  of 
knowledge  of  excitation  mechanisms  and  rates  for  excitation  of  the  state  of  NO.  The  group  of 
7  mechanisms  discussed  here  was  chosen  from  those  reports,  and  the  initial  set  of  rates  that  was 
compared  with  the  present  experimental  data  was  taken  from  Reference  1 1 .  These  rates  are  written 
in  the  form 

kf  =  ATBe-Cn-  (27) 

where  the  values  of  A.  B.  and  C  are  given  in  Table  3.  Ihe  units  of  the  rate  constants  are 
(cm^/molecule  sec)  or  (cm^/molecule*sec). 


A 

B 

C 

1. 

NO-^e 

— » 

N0(A)4€ 

1.2  X  10-10 

0 

63,510 

2. 

N-i-Oh-M 

NO(A)-»-M 

2.47  X  10-31 

-1.24 

0 

3. 

N2-»-e 

N2(A)-Ke 

1.98  X  10-9 

0 

71,586 

4. 

N2-^N 

N2(A)-(-N 

1.53  X  10-2 

-2.23 

71,586 

5. 

N-t-N-i-M 

N2(A)-»-M 

1.27  X  10-27 

-1.60 

0 

6. 

NO(X)-kN2{A) 

NO(A)-^N2(X) 

inf. 

— 

— 

7. 

NCKA) 

NO(X)-Hhv 

5x  106 

— - 

— 

TABLE  3.  INITIAL  RATES  USED.  REFERENCE  1 1 
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The  rates  for  the  reverse  reactions  are  taken  to  be  related  through  the  equilibrium  constants. 
Reaction  7  proceeds  in  only  one  direction  because  the  gas  is  optically  thin. 

An  important  conclusion  of  the  study  of  the  equilibrium  UV  radiation  was  that  the  radiation 
was  not  collision  limited.  This  means  that  at  equilibrium  the  term  k7  (Equation  (21)  Section  3.4)  is 
negligible  compared  with  the  other  terms  in  that  equation.  Typical  values  of  these  other  terms  at 
equilibrium,  for  various  shock  conditions,  are  shown  in  Table  4.  The  initial  pressure  in  the  shock 
tube  was  2.25  ton  for  all  cases. 

Some  of  the  terms,  in  each  case,  should  be  considerably  greater  that  SxlO^sec-l  in  order  to  avoid 
collision  limiting.  Although  this  is  satisfied  for  the  5%  oxygen  experiments  (the  highest 
temperatures),  it  is  not  satisfied  at  the  lower  temperatures  associated  with  the  22.3%  and  40% 
oxygen  equilibrium  conditions.  Thus  if  the  excitation  process  is  to  be  represented  by  these 
relations,  appreciably  faster  rates  must  be  used. 


The  first  two  rate  columns  in  Table  4  show  that  the  electron  reactions  are  the  slowest  of  all 
the  rates.  This  is  related,  of  course,  to  the  small  number  of  electrons  that  are  present  at  these 
relatively  slow  shock  speeds.  The  temperature  dependence  of  the  electron  concentrations  is  such 
that  these  rates  become  more  important  at  higher  temperatures,  but  it  is  unlikely  that  they  are 
involved  in  the  excitation  being  discussed  here. 

The  last  three  rates  have  comparable  magnitudes,  but  the  direct-excitation  rate,  k4f,  has  the 
strongest  temperature  dependence.  The  value  of  the  last  two  terms,  the  recombination  terms,  are 
very  dependent  on  the  degree  of  equilibrium,  since  they  depend  on  the  square  of  the  atom 
concentrations.  At  conditions  near  the  shock,  where  few  atoms  are  present,  they  serve  to  quench 
the  radiation  by  dissociating  the  excited  state  of  NO  orof  Nj. 

Under  the  assumption  that  changes  in  the  ground-state  chemistry  rates  are  slow  compared 
with  the  electronic  excitation  mechanisms.  Equation  (26)  can  be  used  to  describe  the  population  of 
NO  in  the  A  state,  and  thus  (with  Equation  (6))  can  be  used  to  compare  with  the  signal  received 
from  the  Gamma  radiation.  Then 


0 


watts/cm  2sr 


(28) 
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(29) 


Snoa  -  (Z  {lk4fnN2nN+k2fnNnonM+k5fnMnN  I)/ 

{(k4(nN2nN/pi]  +  Ik2fnNnonM/p2]  +k5fiiMnV/P3l} 


The  electron  excitation  terms,  ki  and  k3,  and  the  emission  term  k7  have  been  dropped  from  the 
relation  as  being  negligibly  small.  Only  the  ratios  of  the  remaining  rate  constants  are  relevant  in 
determining  Snqa.  assuming  that  they  are  large  enough  to  make  k7  negligible.  Thus  there  are  two 
variables  in  fitting  the  calculation  with  the  experiments.  Many  calculations  of  Snoa  were 
performed,  using  the  Camac  ground-state  rates  given  in  Table  1  and  varying  coefficients  Aj  for  the 
k2,  k4,  and  ks  electronic  excitation  rates  of  Table  3  in  an  effort  to  match  the  experimental  profiles. 
No  reasonable  match  could  be  obtained. 

The  calculated  ground-state  chemistry  shown  in  Figure  13  is  reasonably  consistent  with  the 
infrared  radiation  data,  but  employ  unrealistic  ground-state  reaction  rates.  If  this  ground-state 
chemistry  is  used  as  the  basis  for  calculation  of  the  electronic  state  populations,  a  fair  agreement 
there  can  also  be  obtained.  A  set  of  calculated  nonequilibrium  emission  profiles  for  the  matrix  of 
230  nm  runs  is  shown  in  Figure  (14).  (The  experimental  results  are  shown  in  larger  scale  in 
Appendix  E.)  These  calculations  used  a  simplified  table  of  coefficients  for  reactions  of  Table  3,  as 
shown  in  Table  5. 


RATE 

A 

B 

c 

1 

1.2  X  10-10 

0 

63,510 

2 

4.0  X  10-34 

0 

0 

3 

2.0  X  10-9 

0 

71,586 

4 

4.6  X  10-8 

0 

71,586 

5 

5.2  X  10-32 

0 

0 

6 

Infinity 

— 

— - 

7 

5  X  106 

TABLE  5.  COEFFICIENTS  USED  IN  CALCULATION  OF  FIGURE  14 

The  pre-exponential  temperature  dependence  was  omitted  in  all  cases.  The  ratios  of  the  Aj  (e.g. 
A4/A2)  are  determined  by  the  shape  of  the  curves.  The  magnitudes  of  the  coefficients  were  chosen 
so  that  the  rates  would  be  large  compared  with  the  radiative  emission  rate.  Rates  I  and  3  do  not 
affect  the  results. 
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INSB  DETECTOR.  3  MM  DIAM.  CHIP 


.  -I  ^ —  COLD  FILTER.  NOM.  5-5.5  ^M 


APERTURE.  1.2  X  3.6  MM 


CAF2  LENS 


3  MM  DIAM.  IMAGE  OF  DETECTOR  CHIP 


SHOCK  TUBE 


ciriiDP  1  MA.IOR  ELEMENTS  OF  THE  IR  RADIOMETER-SHOCK  TUBE  SYSTEM. 
THE  bLm  ^N  the  SHOCK  tube  IS  WEDGE-SHAPED.  OPTIMIZED  FOR  AXIAL 
RFRDLUTION  AT  3MM  AND  FOR  THROUGHPUT  BY  A  9-MM  WIDTH  ACROSS  THE 
TUBE  AT  THE  LOWER  WINDOW.  MAGNIFICATION  IS  UNITY  FOR  THE  DETECTOR 
CHIP  AND  2.5  FOR  THE  APERTURE. 
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FIGURE  2.  TYPICAL  RAW  DATA  TRACE  OF  THE  IR  RADIOMETER,  ILLUSTRATING 
IMPROVEMENT  IN  SCATTERED  LIGHT  REDUCTION.  SEE  TEXT  FOR  DETAILS. 
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FIGURE  6.  RADIATION  MEASUREMENTS,  230  NM 
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PROFILES.  THE  CALCULATED  PROFILES  (DASHED  LINES)  USE  THE 
ELECTRONIC  EXCITATION  RATES  GIVEN  IN  TABLE  5  AND  THE  CALCULATED 
GROUND  STATE  PROFILES  SHOWN  IN  FIGURE  13.  FULL-SCALE  DATA  RECORDS 
ARE  GIVEN  IN  APPENDIX  E. 
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Abstract 

An  experimental  study  is  described  whose 
objectives  include  the  quantitative  measurement  of  the 
ultraviolet  radiation  behind  strong  shock  waves  in  air, 
the  identification  of  the  radiating  species  and  the 
determination  of  the  mechanisms  and  rates  that  govern 
the  emission.  A  shock  tube  was  used  to  generate  3-4  km/s 
shock  waves  through  air  and  O2-N2  mixtures  at  initial 
pressures  between  0.36-2.25  torr.  The  temporal  gas- 
radiance  profiles  recorded  radiometrically  exhibited 
strong  non-equilibrium  overshoots,  followed  by  decay  to 
steady-state  equilibrium.  The  dependence  of  the  peak 
overshoot  values  on  wave  speed  and  pressure  was 
measured  and  shown  to  scale  with  density.  Emission 
spectra  of  the  gas  both  in  the  overshoot  region  and  the 
equilibrium  region  showed  the  NO  gamma  band  system 
between  210  and  340  nm  to  be  the  dominant  radiation 
source.  The  role  of  Oj  Schumann-Runge  radiation  in 
this  region  w'as  also  studied;  in  pure  no  overshoot 
was  recorded.  Profiles  of  the  measured  air  radiation 
have  been  compared  with  computed  profiles  based  on 
translational  and  vibrational  temperatures,  and  are  shown 
to  be  bracketted  by  these.  Infrared  vibronic  emission 
profiles  were  also  measured  and  shown  to  be  a  description 
of  the  NO  concentration  growth.  Further  kinetic  studies 
are  presently  ongoing. 

1.  Introduction 

The  motivation  for  this  research  centers  upon  the 
capability  to  predict  the  ultraviolet  radiation  from  the 
bow  shock  of  an  ascending  flight  vehicle  at  velocities 
between  3-4  km/s  at  altitudes  between  40-60  km.  It  is 
well  known  that  the  strongly  heated  air  behind  the  bow 
shock  of  a  body  moving  at  hypersonic  velocity  can  give 
rise  to  high  radiative  heating  rates  and  to  significant 
radiative  signatures.  These  effects  are  particularly 
important  to  various  aspects  of  reentry  aerodynamics, 
and  much  of  present  understanding  is  a  result  of  intense 
theoretical  and  experimental  efforts  in  this  area  that 
were  undertaken  during  the  late  50's  and  early  60's.  The 
significance  of  each  effect  varies  with  altitude,  velocity 
and  body  configuration,  and  the  capability  to  predict 
these  effects  by  computational  modeling  has  been  an 
ongoing  challenge,  both  in  terms  of  the  physical-chemical 
model  and  in  the  availability  of  data  for  input  parameters 
and  for  comparison  with  predicitons.  The  major  objective 
of  this  study  was  to  provide  experimental  data  that 
directly  supports  predictive  computational  modeling 
efforts  relating  to  non-equilibrium,  bow-shock  radiation. 

The  non-equilibrium  aspect  of  the  radiative 
signature  arises  from  the  fact  that  at  higher  altitudes 
the  decreased  density  and  collision  frequency  of 
molecules  behind  the  shock  front  slows  both  the  chemical 
reactions  that  establish  species  concentrations  and  those 
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that  govern  equipartition  among  internal  energy  states: 
electronic,  vibrational  and  rotational.  Thus,  under  some 
conditions  there  exist  strong  gradients  behind  the  shock 
front,  consisting  of  species  concentrations  and  excited 
state  populations  far  from  equilibrium  with  the  local  (and 
changing)  translational  temperature.  To  describe  this 
non-equilibrium  and  to  predict  its  dependence  on  velocity 
and  altitude  requires  that  the  relevant  kinetic  rates  and 
mechanisms  be  known.  It  is  to  those  that  relate  to 
ultraviolet  and  visible  radiation  that  this  investigation  is 
primarily  directed. 

In  the  more  recent  past,  consideration  of  vehicles 
such  as  the  aeroassisted  and  aerobraked  Orbital  T ransfer 
Vehicles  led  to  computer  code  developments  that  address 
flow  field  predictions  at  higher  altitudes  and  velocities 
(10  km/s)*"^.  Their  application  to  the  lower  velocities 
and  altitudes  of  interest  to  the  present  problem  is  not 
clear,  owing  to  the  widely  different  kinetic  chemistry 
involved.  The  need  for  an  experimental  data  base  over 
the  relevant  velocity-altitude  range  constitutes  the  major 
objective  of  these  studies. 

The  most  noteworthy  and  interesting  difference 
between  early  shock  tube  studies  and  those  reported 
herein  lies  in  the  newer  instrumentation  available  today. 
Most  especially  useful  were  a  multi-channel,  high-speed 
digital  data  recording  system,  and  an  optical  multichannel 
analyzer.  This  electronic  spectrometer  permitted  short, 
precisely  synchronized  exposures  to  be  made  with  good 
wavelength  resolution. 

The  shock  tube  and  the  instrumentation  is  briefly 
described  below.  Typical  data  and  their  calibration  are 
presented,  together  with  a  subsequent  discussion  of 
results. 

2.  Experimental 

Shock  Tube 

The  shock  tube  employed  in  this  research  was  used 
in  a  number  of  radiation  and  chemical  kinetic  studies 
over  several  decades.***^ 

A  schematic  diagram  of  the  overall  experiment  is 
shown  in  Figure  1.  The  steel  driven  section  is  30  feet 
long,  3-inch  inside  diameter.  The  driver  section  is  5 
feet  long,  and  was  originally  designed  for  operation  at 
15,000  psi  using  pure  hydrogen  as  the  driver  gas.  To 
date  on  this  program  a  maximum  of  4,000  psi  was  used; 
most  tests  only  required  1,000  psi  of  H2,  achieved  by 
means  of  press-scribed,  0.035-inch-thick  steel 
diaphgrams.  The  routine  double-diaphragm  technique* 
provided  excellent  run-to-run  reproducibility  in  wave 
speed  (~1“6>  and  radiation  records. 
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permitted  spectra  of  both  equilibrium  and  nonequilibrium 
air  radiation  to  be  obtained.  This  instrument  was  not 
calibrated  absolutely,  although  the  overall  responsivity 
as  a  function  of  wavelength  was  measured  and  taken  into 
account  for  each  grating. 

A  HgCdTe  detector-filter  radiometer  was  also 
deployed  to  measure  the  MO  fundamental  vibronic  band 
at  5  pm,  to  obtain  data  relevant  to  the  NO  vibrational 
kinetics.  The  infrared  signa!-to-noi$e  ratios  were  low, 
but  adequate  radiance  profiles  were  obtained  and  will  be 
discussed  later. 

An  uncalibrated,  wider  bandpass,  1  mm-wide-beam 
radiometer  at  250  nm  provided  redundant,  primarily 
housekeeping,  data. 


Fig.  1  Principal  elements  of  the  experiment. 

Vacuum  levels  achieved  using  diffusion  pumps  with 
cryogenically  cooled  traps  were  routinely  in  the  10'^  torr 
range  with  negligible  leak/outgassing  rates  over  the  15 
minute  pre-test  gas  loading  interval.  Initial  test  gas 
pressures  were  measured  to  better  than  l<?i  with  an  MKS 
Baratron  unit.  Scientific  grade  O2  and  N2 

were  used  throughout,  with  H2O  and  CO2  at  less  than 
1  ppm.  After  each  test  the  tube  was  scrubbed  with 
alcohol. 

Thin-film  heat  transfer  gauges  were  used  to  record 
shock  arrival  at  seven  sections  along  the  tube.  These 
signals  were  differentiated  and  combined  for  digital 
recording.  This  record  then  provided  a  common  time 
line  for  all  subsequent  measurements,  as  discussed  in  the 
following  section. 

Instrumentation 


OHC  CO  Tt  DCTECTOR 

" 


An  instrumentation  station  with  four  ports  was 
located  about  six  feet  from,  the  end  of  the  shock  tube. 
Figure  7  illustrates  the  four  measurements  that  were 
made  at  that  station  for  the  test  results  reported  in  this 
paper.  The  primary  quantitative  measurements  were 
made  using  a  radiometer  comprising  a  UV-quartz  window, 
a  IP28  photomultiplier  tube,  an  interference  filter  (260 
nm,  10  nm  FW'HM,  I5‘*('  transmission)  and  two  slits  to 
define  the  optical  beam  to  3  mm  maximum  width.  These 
optical  parameters  were  selected  after  a  number  of 
exploratory  tests  were  made  to  establish  spectral  purity, 
and  to  assess  signal  strengths  and  spatial  and  temporal 
radiance  gradients.  This  radiometer  was  calibrated  by 
using  an  Optronics  standard  irradiance  deterium  arc  lamp, 
MBS  traceable.  In  the  calibration  setup,  a  rotating  mirror 
of  known  reflectance  provided  a  dynamic  lamp  scan  using 
the  same  time  constant  and  slit-PMT  cathode 
configuration  as  in  the  experiment.  Calibrations  slightly 
overlapping  the  range  of  signal  levels  of  the  experiment 
were  made,  but  were  limited  by  low  lamp  signals  through 
the  narrow  bandpass  filter.  Linearity  of  output  over  the 
entire  range  was  confirmed  by  using  the  detector 
unfiltered  in  an  inverse-square/distance  set  of 
measurements. 


Fig.  2  Instrumentation  at  the  measurement  station. 


Finally,  it  should  be  noted  that  another  radiometer 
was  deployed  at  a  more  downstream  station  in  the  shock 
tube,  its  filter  at  307  nm  was  selected  to  monitor  the 
10,  0)  transition  of  the  OH  (A-Xt  band  system.  Its  purpose 
was  to  monitor  shock  tube  performance  by  recording  the 
passage  of  the  reacting  H2-driver/air  driven  gas 
interface.  It  thereby  confirms  the  length  (or  duration) 
of  the  shock-processed  test  gas  "slug”,  which  ranged  from 
60-100  microseconds. 

The  OMA  was  a  completely  self-contained  system; 
a  computer  was  used  for  synchronization,  exposure 
control  timing,  data  recording  and  processing.  All 
radiometer  signals  and  the  wavespeed  data  were  fed  to 
a  datalab,  8-channel  digital  recorder  with  a  2  MHz 
sampling  rate.  This  arrangement  permitted  all  shock 
tube  test  data  to  be  correlated  to  within  0.5  microsecond 
resolution;  they  included  the  radiometer  outputs  at  two 
gain  settings,  the  wavespeed  time-of-arrival  pulses,  and 
a  signal  from  the  OMA  that  permitted  the  exact  exposure 
time-interval  of  the  spectrometer  to  be  recorded. 


An  optical  multichannel  analyzer  (EG  8cC(PARC 
OMA  3'  was  also  deployed  as  shown.  Its  main  function 
was  to  provide  emission  spectra  throughout  the  UV  and 
visible  wavelength  range,  primarily  to  assess  the  spectro- 
chemica!  purity  in  the  shock  tube  test  gas,  and  to  verify 
the  origin  of  the  radiation  recorded  by  the  radiometers. 
This  instrument  proved  to  be  highly  effective;  its 
operational  feature  of  precise  electronic  shuttering 


3.  Typical  Data 


This  section  will  illustrate  the  nature  of  the 
various  data  derived  from  the  experiment,  the  techniques 
used  for  data  reduction  and  calibration,  and  how  the 
records  have  been  interpreted. 
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The  wave  speed  data  format  comprises  a  series 
of  spikes  that  corresponds  to  shock  arrival  along  about 
half  the  shock  tube  length.  The  intervals  between  spikes 
could  be  measured  to  0.5  microseconds,  from  which  a 
plot  like  Figure  3  was  generated  for  each  run.  A  second- 
order  curve  was  then  fitted  to  the  data,  and  its  slope 
at  the  measurement  station  was  taken  as  the  shock  wave 
velocity  for  the  test.  The  inputs  for  all  shocked-gas 
property  calculations  were  the  initial  test-gas  pressure, 
the  shock  velocity  and  the  ambient  temperature. 


Fig  3  Reduced  wave  speed  data  Measurement 
station  at  2843  mm. 

In  Figure  it  are  shown  very  careful  tracings  of  the 
raw-data  radiance  profiles  obtained  with  the  2fC  nm 
radiometer  for  two  tests  in  air  at  the  cited  conditions. 
For  all  conditions  of  interest,  the  overall  profiles  were 
similar,  namely,  a  strong,  rapidly  rising  overshoot  region, 
followed  by  a  decay  to  a  steady-state  equilibrium.  The 
features  of  interest  are  the  time-to-peak,  peak  signal 
levels,  the  decay  rate,  and  the  signal  levels  at 
equilibrium. 

As  will  be  shown  below,  this  overshoot  is 
dominated  by  200-320  nm  emission  from  the  gamma  band 
system  of  NO.  The  NO  begins  to  form  at  the  onset  of 
O2  dissociation  behind  the  shock  front.  Rapid  NO 
electronic  excitation  takes  place,  owing  to  the  high 
translational  temperature  immediately  behind  the  shock 
front.  This  temperature  decays  rapidly,  as  it  is  primarily 
dominated  by  the  O2  dissociation  rate.  Thus,  the 
electronic  excited  state  population  is  determined  by  a 
number  of  competing  and  coupled  kinetic  processes, 
whose  result  is  the  measured  observables  in  Figure  4. 
The  chemical  kinetic  mechanisms  and  rates  involved 
constitute  a  major  objective  of  this  research. 

One  approach  to  this  objective  is  to  measure  the 
overshoot  profiles  over  a  broad  range  of  N2/O2  mixture 
ratios,  pressures  and  temperatures.  One  example  is 
presented  in  Figure  5,  in  which  pure  O2  was  used  as  the 
test  gas.  In  this  case  the  radiation  was  shown  to  arise 
from  the  Schumann-Runge  band  system  of  O2.  The 
striking  difference  in  this  profile  is  the  lack  of  any 
overshoot.  This  feature  of  O2  understood^.  It 

results  from  a  dominating  radiationless  depopulation  path 
of  the  excited  state  to  atoms;  and  the  atom  density- 
increases  monotomically,  without  an  overshoot.  This 
feature  is  significant  in  the  kinetics  of  the  overshoot. 
If  the  same  kinetics  prevail  in  air,  then  the  role  of 
O2  in  the  overshoot  will  be  limited  by  its  value  in  the 
equilibrium  region  of  the  profile  of  Figure  4. 


Typical  5  pm  infrared  profiles  are  shown  in 
Figure  6.  Three  profiles  were  superimposed;  all  are  at 
initial  air  pressures  of  2.25  torr,  with  shock  velocities 
at  3.22,  3.50  and  3.83  mm/ps.  Despite  the  much  lower 
SNR,  the  comparisons  show  a  definite  trend  of  overshoot 
behavior.  The  interpretation  and  role  of  the  overshoot 
is  addressed  in  Section  4. 

The  DMA  has  proven  to  be  invaluable  to  this  study 
by  providing  emission  spectra  of  the  shock-heated  gases. 
Initially,  broad  wavelength  coverage  spectra  were  taken 
to  identify  the  dominant  radiating  species/bands,  and  to 
assess  overall  system  purity  levels. 


torr,  3.50  mm/ps. 
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Fig.  6  Selected  Sum  radiation  profiles  illustrating  the 

dependence  on  shock  velocity.  Initial  air  pressures: 
2.25  torr;  velocities  in  mm/us:  (a)  3.83,  (b)  3.50. 
(c)  3.22. 


Primary  concern  of  spectral  priority  is  the 
verification  that  the  quantitative  signals  measured  by 
the  radiometers  are  indeed  fromthe  species  of  interest. 
A  typical  low-resolution  spectrum  of  shock-heated  air  is 
shown  in  Figure  7.  In  this  case  the  spectrum  was 
integrated  over  30  us  during  the  equilibrium  portion  of 
a  profile  such  as  in  Figure  <i(a),  between  50-80  ps.  For 
the  profile  plots,  time  t  s  0  is  measured  from  a  time- 
of-arriva!  gauge  located  73  mm  ahead  of  the  a-port 
instrumentation  section  centerline. 

The  recorded  spectrum  clearly  shows  the  band 
structure  of  the  NO  gamma  band  system  in  both  the 
first  and  second  order  of  the  grating.  No  hydrocarbon 
impurity  bands  are  in  evidence.  To  date,  only  the 
ubiquitous  Na  D  line  and  the  306.<i  nm  band  of  the  OH 
(0,  01  A-X  have  been  seen.  The  origin  of  the  OH  has 
not  been  pursued;  it  is  of  variable  signal  strength,  and 
felt  to  arise  from  HjO  on  the  walls  of  the  tube.  We 
note  that  with  hydrogen  driven  gas,  combustion  at  the 
interface  occurs  in  every  test.  Fortunately,  unambiguous, 
quantitative  radiometer  measurements  at  wavelengths 
below  280  nm  can  be  made  despite  the  trace  OH  presence. 
To  date,  its  role  in  kinetics  or  energetics  has  been 
ignored. 

The  presence  of  On  radiation  in  the  spectrum  of 
Figure  7  is  clearly  overwhelmed  by  the  NO  system.  A 
spectrum  of  shock-heated  O2  at  higher  resolution  is  shown 
in  Figure  8.  As  can  be  seen,  the  spectrum  below  30D 
nm  is  fairly  featureless,  being  made  up  of  many 
overlapping  bands.  Identification  of  several  most  obvious 
bands  is  given  at  longer  wavelengths.  The  spectra  of 
Figures  7  and  8  can  be  compared  with  the  computed 
spectra  of  Figures  12  and  13. 


Fig.  7  Typical  low-resolution  spectrum  of  shock-heated 
air  in  equilibrium  at  3800<2K. 
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Fig,  8  High  resolution  spectrum  of  oxygen  in  equilibrium 
at  3380OK. 


It.  Results 


Nonequilibrium  Radiation 

Velocity-Altitude  Scaling 

The  nonequilibrium  radiation  from  the  shock- 
heated  gas  was  measured  for  a  number  of  shock  velocities 
at  each  of  three  initial  shock-tube  pressure.  The  initial 
air  pressures  of  2.25,  0.66  and  0.36  torr  were  chosen  to 
correspond  to  flight  altitudes  of  40,  50,  and  55  km. 
Lower  pressures  were  avoided  on  the  usual  grounds  of 
difficulties  with  boundary  layers,  shock  curvature,  etc. 

When  the  correspondence  between  the 
nonequilibrium  zone  behind  the  shock  wave  in  the  tube 
and  that  behind  a  blunt-body  bow  shock  is  examined,^*^^ 
it  can  be  seen  that  the  peak  overshoot  value  is  the  most 
critical  factor  in  determining  the  radiation  from  the  gas 
cap.  The  region  of  the  decay  to  equilibrium  and  the 
equilibrium  radiance  arc  compressed,  lie  close  to  the 
body,  and  are  contained  in  the  boundary  layer. 
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The  results  of  some  20  runs  are  plotted  in  Figure 
9,  in  which  the  peak  overshoot  radiance  values  were 
analyzed  by  the  densities  at  the  altitudes  shown.  The 
direct  scaling  with  density  is  readily  apparent,  as  is  the 
factor  of  100  in  intensity  over  the  velocity  range  shown. 
To  obtain  the  absolute  radiance  units,  the  signals  from 
the  radiometer  were  directly  ratioed  to  those  from  the 
known  Dj  standard  lamp,  with  account  taken  of  the  shock 
tube  diameter  and  the  optica)  throughput  (area  x 
steradiance)  of  the  radiometer.  It  should  be  mentioned 
that  the  measurements  shown  are  in-band,  corresponding 
to  the  260  nm  filter.  In  order  to  extrapolate  these  data 
to  other  wavelength  intervals,  the  relevant  spectra  of 
NO  and  Oj  at  the  conditions  of  interest  would  need  to 
be  considered.  For  this  purpose  the  spectra  derived  with 
the  OMA  were  particularly  significant. 


VELOCITY  (KW'S) 

Fig.  9  Velocity  and  altitude  dependence  of  the  peak 
overshoot  value  of  the  in-band  volumetric  air 
radiance.  Data  normalized  by  density  ratios. 


Electronic  State  Kinetics 


The  ultimate  goal  of  the  kinetics  studies,  both 
theoretical  and  experimental,  is  to  determine  the 
mechanisms  and  rates  for  the  population  of  the  excited 
NOr^'  state.  Except  for  the  case  of  NOfA'  formation 
directly  from  N  and  O  atoms,  the  most  likely  other 
mechanisms  involve  collisions  with  the  NO  molecules 
being  formed  behind  the  shockfront.  The  computational 
aporoach  taken  to  data  is  to  use  a  kinetic  code  to 


calculate  the  formation  of  the  ground  state  chemical 
species  first.  This  implies  that  excited  electronic  state 
kinetics  can  be  decoupled.  This  convenience  then  permits 
separate  interations  on  these  mechanisms  and  rates, 
without  recalculation  of  all  species  concentrations. 

The  first  step  has  been  achieved  by  making  use 
of  an  updated  version  of  Calspan's  Kinetic  Normal  Shock 
(KNS)  code,!'  redesigned  to  run  on  an  IBM  PC  AT.  It 
computes  the  thermal  and  chemical  relaxation  of 
mixtures  of  reacting  gases  behind  normal  shock  waves. 
The  menu-driven  program  enables  the  user  to  produce  a 
convenient  listing  of  selected  parameters  versus  particle 
time  and  distance  on  a  file  for  printing  and  generate 
plot  files  which  may  be  viewed  on  the  screen  and/or 
routed  to  a  printer. 

Table  1  presents  a  listing  of  the  II  reactions  and 
the  rates  used  for  the  computations  for  air. 

Table  1 

REACTION  RATE  COEFFICIENTS 
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An  example  of  some  results  obtained  with  these 
reaction  rates  is  shown  in  Figure  10.  Separately  shown 
are  the  translational  and  N2  vibrational  temperatures, 
and  the  NO  and  O2  concentrations  behind  the  shock  front. 
In  Figure  1 1  the  computed  excited  fA>  state  population 
of  NO  is  plotted,  assuming  local  equilibrium  with  each 
of  the  temperatures.  For  clarity  in  comparing  the 
temporal  features  of  the  overshoots,  the  data  were 
normalized  to  unity  peak  values.  Also  shown  is  the 
normalized  radiometer  record  for  the  same  conditions. 
It  can  be  seen  that  the  excitation  rate  is  characterized 
by  a  temperature  intermediate  to  those  shown. 

Further  efforts  are  presently  underway  to  obtain 
similar  data  over  a  broad  range  of  O2/N2  mixture  ratios, 
and  to  use  this  extended  data  base  in  formulating  the 
operant  mechanisms  and  rates  that  determine  the 
radiative  overshoot  profiles. 

Infrared  Vibrational  Kinetics 


The  infrared  radiation  of  NO  in  air  has  been 
measured  for  a  series  of  shock  speeds  and  for  an  initial 
shock-tube  pressure  of  2.25  torr.  The  emission  results 
are  shown  in  Figure  6  for  three  experiments  spanning 
this  range.  The  spatial  resolution  of  the  optical  system 
is  about  5  mm,  or  on  the  scale  shown  in  Figure  6,  about 
l.ii  usee.  Vibrational  relaxation  times  for  NO  fin  Argon' 
range  from  30  to  fcO  usec'^>'^  (particle  time'  under  these 
conditions  of  temperature  and  pressure,  and  are  probably 
considerably  shorter  in  air  (they  are  50  times  shorter  in 
pore  NO'^>*^).  The  laboratory  time  Ithe  parameter  used 
in  Figure  6'  is  shorter  than  the  particle  time  by  a  factor 
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Fig.  10  Direct  output  plots  of  the  computed  temperature 
and  species  concentration  gradients  behind  the 
shock  front.  Velocity;  3  84  Km  s;  initial  air  pressure: 
2.25  torr. 


Fig.  11  Normalized  NO  (A)  state  population  profiles 
in  local  equilibrium  at  (a)  translational  tem¬ 
perature  (b)  N  ^  vibrational  temperature 
(c)  260  nm  radiometer  measured  profile. 

Same  test  conditions  as  for  Fig.  10. 

equal  to  the  density  ratio  across  the  shock,  which  is 
between  6  and  10  for  these  experiments.  Thus,  the 
vibrational  relaxation  time  of  NO  is  of  the  order  of 
several  microseconds  in  laboratory  time.  The  observed 
intensity  is  then  proportional  to  the  an  7unt  of  NO  present 
and  to  a  weak  function  of  the  vibrational  temperature 
of  the  NO.  For  practical  purposes  the  position  of  the 
intensity  maximum  would  be  expected  to  occur  close  to 
the  maximum  value  of  the  NO  concentration.  The 
measured  times  to  maximum  are  in  fact  longer  than 
those  calculated  with  the  standard  11-reaction  model 
given  in  Table  1.  The  calculated  NO  maxima  occur  at 
1.9  iisec,  7.S  ijsec  and  36.5  usee  for  the  experiments 
shown  in  Figure  6,  compared  with  experimental  values 
of  about  S  usee,  22  usee  and  as  usee.  One  explaination 
would  be  a  slower-than-expected  relaxation  of  NO 
vibration.  The  slow  decay  of  NO  concentration  at  times 
after  the  maximum  is  reached  would  account  for  the 
final  decay  of  the  radiation. 

Equilibrium  Radiation 


NO  and  O'y  f  Values 

The  time-dependent  measurements  of  the  radiation 
at  260  nm,  as  shown  in  Figure  4,  can  be  used  to  obtain 
an  equilibrium  radiation  value  at  long  times,  and  these 
can  be  compared  with  previous  equilibrium  experiments 
by  evaluation  of  the  optical  f  values  of  the  radiating 
systems.  The  experimental  data  were  reduced  by  using 
a  numerical  calculation  which  approximates  the  spectral 
distribution  by  replacing  the  actual  vibrational-rotational 
structure  with  a  simple  smoothed  Q-branch  structure  fsee 
Figures  12  and  13)'**,  and  performing  a  numerical  integral 
of  the  intensity  multipled  by  the  filter  function.  If  the 
calculation  is  performed  for  O2  and  NO,  using  f  values 
of  unity  and  number  densities  of  unity  for  each  species, 
reference  values  for  the  in-band  radiation,  and  I°\Ot 
are  obtained.  The  in-band  radiation  to  be  expected  from 
air  would  then  be 

Tatr  *  ^2  fo3  I02  *’no  fNO  Ino 


or 


ho2  ^02 


'03 


'NO 


Since  1°no  '°C2  known  from  the  calculation, 
and  n>v;(>  and  n^j  are  knov-  from  the  equilibrium  shock 
calculation,  the  experimental  value  of  Igj^  divided  by 
n^2  '°02  plotted  vs  niv;(^lOnj(^/nf>2l°C'2- 

intercept  and  slope  of  the  straight  line  yield  the  i-values 
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of  and  NO.  The  values  obtained  from  eleven 
experiments  in  O,  and  air  are  fQ2  =  0.052  and  z 
0.0091  at  a  wavelength  of  260  nm,  as  shown  in  Figure 

14.  The  value  for  Oj  is  somewhat  lower  than  that  sug¬ 
gested  by  Keck  et  al*^  (0.075)  or  by  Krindach  et  a!'^ 
10.091  for  this  wavelength.  The  NO  value  is  consistent 
with  the  values  given  by  Schadee^^. 


wavelength  (nm) 

Fig  12  Nitric  oxide  gamma  band  system  (4000K) 


References 

1.  Park,  C.,  Calculation  of  Nonequilibrium  Radiation 
in  the  Flight  Regimes  of  Aeronautical  Orbital 
Transfer  Vehicles,  Paper  84-0306  at  AlAA  22nd 
Aerospace  'l-iences  Mtg.,  Reno,  3an.  1984. 

2.  Park,  C.,  Problems  of  Rate  Chemistry  in  the  Flight 
Regimes  of  Aeroassisted  Orbital  T ransfer  Vehicles, 
Paper  84-1730  at  AlAA  I9th  Thermophysics 
Conference,  Snowmass  Colo.,  June  1984. 

3.  Candler,  G.  and  MacCormack,  R.,  The  Computa¬ 
tions  of  Hypersonic  Ionized  Flows  in  Chemical  and 
Thermal  Nonequilibrium  Paper  88-0511  at  AlAA 
Aerospace  Sciences  Mtg.,  Reno,  Jan.  1988. 

4.  Wurster,  W.H.,  Measured  Transition  Probability  for 
the  First-Positive  Band  System  of  Nitrogen,  J. 
Chem.  Phys.  36,  April  1962. 

5.  Marrone,  P,.V.  and  Wurster,  W.H.,  Measurement 
of  Atomic  N  and  C  Using  Shock  Tube  VUV 
Spectrometry,  JQSRT  21  327  1971, 

6.  Rothe,  D.E.,  Radiative  lon-EIectron  Recombination 
J.  Quant.  Spectres.  JQSRT  9,  49,  1969. 


1.00 


L.J  i.n,  M  ■■  I  i  ^  . 

I _ _  : 

L.  .  1  1  I  i  ' 

□r  ' 

iLLJ _ -  i  . 

,  :  .Li,  1  f-  j 

'  i 

_ .  '  ' 

'  :  i  " 

1 

!  1 

t  t  1 

1  t  1  1 

L_! _ ^ _ L_ 

'III 

; 

200  240  280  320  360  400 


wavelength  (nm) 

Fig  i3  Oxygen  Schumann-r  ,nge  band  system  {3383K) 


Fig  14  Intensify  measurements  in  oxygen  and  in  air. 


7.  Boyer,  D.V',,  Rate  Coefficient  Determinations  of 
Important  Reactions  Involving  HBO2  and  OBF, 
AFRPL  TR-S3-045,  March  1984. 

8.  Herzberg,  C.,  Spectra  of  Diatomic  Molecules,  D. 
Van  Nostrand  Co.,  Princeton,  N.J.,  1955. 

9.  Cibson,  Tt'.E.  and  Marrone,  P.V.,  "A  Correspond¬ 
ence  Between  Normal  Shock  and  Blunt-Body  Flows" 
CAL  QM-1626-A-7,  June  1962. 

10.  Park,  C,,  Assessment  of  a  Two  Temperature 
Kinetic  Model  for  Ionizing  Air,  AlAA  Paper  87- 
1574,  Hawaii,  1987, 

1 1.  Marrone,  P.V.,  et  al,,  Inviscid,  Nonequilibrium  Flow 
Behind  Bow  and  Norma!  Shock  Waves,  Cornell 
Aeronautical  Laboratory  Report  No.  QM-1626-A- 
12  (1-111)  May  1963,  October  1966. 

12.  Wray,  Kurt  L.  A  Shock  Tube  Study  of  the 
Vibrational  Relaxation  of  Nitric  Oxide,  J.  Chem 
Phys  36  2597  (1962). 

13.  Kamimoto,  Goro  and  Matsui,  Hiroyuki.  Vibration 
Relaxation  of  Nitric  Oxide  in  Argon,  J.  Chem  Phys 
53  3987  (1970). 

14.  Treanor,  C.E.  and  W'illiam*  M.J.  A  Method  for 
Calculating  Diatomic  Spectra  Using  a  Digital 
Computer  CAL  Report  QM1626-A-5,  1962. 

15.  Keck,  J.C.,  Allen,  R.A.  and  Taylor,  R.L. 
"Electronic  Transition  Moments  for  Air  Molecules", 
JQSRT  3,  33>  (1963). 

16.  Krindach,  N.I.,  Sobolev,  N.N.,  and  Tunitskii,  L.N. 
Determination  of  the  Electronic  Transition 
Moments  for  the  Schumann-Runge  Bands  of  O2, 
Opt.  Spectrosc.  (Eng.  trans.1  25,  ?26  (1963). 

17.  Schadee,  Aert.  The  Relation  Between  the 
Electronic  Oscillator  Strength  and  the  Wavelength 
for  Diatomic  Molecules,  JQSRT  7  169  (19671. 


A-7 


UNCLASSIFIED 


SIFICATlON 


MASTER  COPY 


FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


2a.  SECURITY  OASSlElCATION  AUTHORITY 


2b.  DECLASSIFICATION /OOWNGRAOING  SCHEOUL 


4.  PERFORMING  ORGANIZATION  REPORT  NUMIER($) 


lb.  RESTRiaiVC  MARKIN 


2  .  DtVTRlPU'ilON /  availability  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


fa.  NAME  OF  PERFORMING  ORGANIZATION 

CALSPAN-UB  Research  Center 


6c  ADDRESS  (CKy,  Smn,  sndZmCodt) 

Box  400 

Buffalo,  N.y.  14225 


fa.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

U.  S.  Army  Research  Office 


Be  ADDRESS  (Oty,  Sum.  and  ZIP  Code) 


6b  OFFICE  SYMBOL 
(ff  sppUubh) 


Bb.  OFFICE  SYMBOL 
pf  afipUesbh) 


7a.  NAME  OF  MONITORING  ORGANIZATION 
U.  S.  Army  Research  Office 


7b.  ADDRESS  (Oly,  State,  and  ZIP  Cede! 

P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


B.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


to  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
CLEMENT  NO 


PROJEa 

NO 


P.  0.  Box  12211 

Rese..rch  Triangle  Park,  NC  27709-2211 


1 1 .  title  (Includt  Stcurity  a*tufi<»tion) 

Nonequilibrium  UV  Radiation  and  Kinetics  behind  shock  waves  in  air 


t2  PERSONAL  AUTHOR(S) 

_  _  W.H.  Wurster,  C.E.  Treanor  and  M.J.  Williams 


1 3a.  TYPE  OF  REPORT  It  3b  TIME  COVERED  |t4.  DATE  OF  REPORT  (Year,  MOntA  Day)  ll  5  PAGE  COUNT 

Reprint  I  FROM _ TO _  I  I 


16  SUPPLEMENTARY  NOTATION  „  ,  ,  , 

The  view,  opinions  and/or  findings  contained  In  this  report  are  those 

as.  an  official  Department  of  the  Army  position, 


17  COSATI  CODES  I  IB.  SUBJECT  TERMS  ^ontinwe  on  revene  H  necesary  and  idtntify  By  Bloc*  numb*r) 

f'ELO  I  GROUP  I  SUBGROUP  I  Radiation  Shock-Heated  Air 

Ultraviolet  Nonequilibrium  Air 


19  ABSTRACT  (Continue  on  /eveiw  if  necesary  and  Uttttify  By  bhdi  numBer) 

An  ejtperimental  study  is  described  whose  objectives  include  the  queintitative  measurement  of 
the  ultraviolet  radiation  behind  strong  shock  waves  in  air,  the  identification  of  the 
radiating  species  and  the  determination  of  the  mechanisms  and  rates  that  govern  the  emission 
■  A  shock  tube  was  used  to  generate  3-4km/s  shock  waves  through  air  and  O2-N2  mixtures  at 
initial  pressures  between  0.36-2.25  torr.  The  temporal  gas-radiance  profiles  recorded 
radiometrically  exhibited  strong  non-equilibrium  overshoots,  followed  by  decay  to  steady- 
state  equilibrium.  The  dependence  of  the  peak  overshoot  values  on  wave  speed  and  pressure 
was  measured  and  shewn  to  scale  with  density.  Bnission  spectra  of  the  gas  both  in  the  over¬ 
shoot  region  eutd  the  equilibrium  region  showed  the  NO  ganiTB  beind  system  between  210  and 
340  nm  to  be  the  dominant  radiation  source.  The  role  of  O2  Schumann-Runge  radiation  in  this 
region  was  also  studied;  in  pure  O2  no  overshoot  was  recorded.  Profiles  of  the  measured  air 
radiation  have  been  compared  with  computed  profiles  based  on  translational  and  vibrational 
temperatures,  eind  are  shewn  to  be  bracketted  by  these.  Infrared  vibronic  emission  profiles 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRAa 
□  UNCLASSIFIEDAJNLIMITED  □  SAME  AS  RPT. 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 


121.  ABSTRAa  SECURITY  CLASSIFICATION 

Unclassified  _ 


22b  TELEPHONE  (BkIinN  A/«a  Cede;  22c.  OFFICE  SYMBOL 


DD  FORM  1473.64  MAR 


B3  APR  tdttion  may  b«  uMd  until  eKhaunad 
Ail  other  editiont  aro  obiolata 


UNCLASSIFIED 


UNCLASSIFIED 

■■CUIMTV  CCAMiriCATlOW  OF  TWH  _  _ 

I  19.  (cont.)  Further  kinetic  studies  are  presently  ongoing. 


UNCLASSIFIED 


MCUMITT  CUAWiriCATlON  OP  THIS  PAOC 


AIAA-90-1666 

Kinetics  of  UV  Production  Behind  Shock 
Waves  in  Air 

W.  H.  Wurster,  C.  E.  Treanor  and  M.  J.  Williams 
Calspan-UB  Research  Center, 

Buffalo,  NY 


AIAA  21st  Fluid  Dynamics,  Plasma  Dynamics 
and  Lasers  Conference 

June  18-20,  1990  /  Seattle,  WA 


r 


For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  and  Astronautics 
370  L'Enfant  Promenade.  S.W.,  Washington,  D.C.  20024 


KINETICS  OF  UV  PRODUCTION  BEHIND  SHOCK  WAVES  IN  AIR* 


W.H.  Wurster*.  C.E.  Treanor**,  and  M.J.  Williams*** 
Calspan-UB  Research  Center 
Buffalo.  New  York 


Abstract 

In  a  series  of  shock-tube  e}q>eriments,  the 
time-dependent  ultraviolet  radiation  behind  the 
shock  wave  has  been  measured  for  various  Os/Nj 
mixtures  and  for  a  range  of  shock  velocities  horn 
3.24  to  3.86  mm/^sec.  From  the  equilibrium  portion 
of  the  experiments,  the  values  of  the  squared  transi¬ 
tion  moment  for  the  O2  Schumarm-Runge  and  the 
NO  Gamma  band  systems  were  measured  at  230  nm 
to  be  iRelsu*  1-3  a.u.  and  iRel^  «  0.21  a.u.  For  the 
nonequilibrium  region  directly  behind  the  sho^k. 
calibrated  time  histories  of  NO  Gamma  radiation  a:  e 
presented.  The  kinetic  mechanisms  governing  this 
nonequilibrium  radiation  have  been  investigated,  but 
values  for  specific  rate  consums  have  not  yet  been 
determined.  Concurrent  infrared  measurements 
which  viewed  the  NO  IR  system  at  1850  cm*^  show 
agreement  with  calculated  NO  equilibrium  values, 
but  did  not  have  the  time  resolution  to  follow  the 
details  of  the  nonequilibrium  emission.  Plans  for  im¬ 
proved  instrumentation  and  further  measurements 
are  described. 

1.  Introduction 

This  paper  describes  the  latest  results  of  an  on¬ 
going  research  program  whose  motivation  centers 
upon  the  capability  to  predict  the  ultraviolet  radia¬ 
tion  from  the  bow  shock  of  a  vehicle  ascending  be¬ 
tween  40-60  km  at  velocities  between  3-4  km/s. 
The  objectives  of  the  program  include  providing 
quantitative  measurements  in  suppon  of  a  sounding 
rocket  flight  exp>eriment,  the  identification  of  the  ra¬ 
diating  species,  and  the  dependence  of  the  radiation 
on  altitude  and  velocity.  In  addition,  basic  dau  are 
being  acquired  in  support  of  modelling  this  non- 
equilibrium  radiation. 

The  non-equilibrium  aspect  of  the  radiative 
signature  arises  from  the  fact  that  at  higher  altitudes 
the  decreased  density  and  collision  frequency  of 
molecules  behind  the  shock  front  slows  both  the 
chemical  reactions  that  esublish  species  concentra¬ 
tions  and  those  that  govern  equipartition  among  in¬ 
ternal  energy  sutes:  electronic,  vibrational  and 
routional.  Thus,  under  some  conditions  there  exist 


strong  gradients  behind  the  shock  front,  consisting  of 
species  concentrations  and  excited  stau  populations 
far  from  equilibrium  with  the  local  (and  ^nging) 
translational  temperature.  To  describe  this  non- 
equilibrium  and  to  predia  hs  dependence  on  veloc¬ 
ity  and  altitude  requires  that  the  relevant  kinetic 
rates  and  mechanisms  be  known.  It  is  to  those  that 
relate  to  ultraviolet  and  visible  radiation  that  this  in¬ 
vestigation  is  primarily  directed. 

A  description  of  the  shock  tube  experiment 
and  early  results  have  been  presented  in  Reference 
1.  They  will  be  very  briefly  discussed  below,  for 
overall  completeness.  New  instrumenution  and  an 
extension  of  the  range  of  measurements  are  next  de¬ 
scribed.  They  provide  an  extensive  data  base  whose 
analysis  is  underway,  with  the  goal  of  providing  a 
model  to  describe  the  electronic  sute  kinetics  that 
govern  the  nonequilibrium  ultraviolet  radiation  from 
shock-heated  air.  The  results  to  date  of  these  ongo¬ 
ing  analyses  are  discussed,  and  present  plans  to  com¬ 
plete  the  program  are  outlined. 

2.  Review  of  Previous  Work 

It  is  especially  convenient  that  the  flight  condi¬ 
tions  of  interest,  namely,  shock  velocities  of  3-4 
km/s  at  altitudes  above  40  km  (2.25  torr)  can  be 
exactly  duplicated  in  a  shock  tube.  Thus,  the  meas¬ 
urement  of  the  radiation  profile  as  the  shock-heated 
air  passes  a  given  station  can  be  readily  and  simulu- 
neously  achieved  for  a  number  of  wavelengths  and 
with  different  instruments.  Further,  it  can  be  shown 
that  a  profile  measured  in  the  shock  ribe  can  be 
transformed  to  that  behind  the  bow  shock  wave.  The 
use  of  this  “correspondence”  then  permits  the  direct 
use  of  the  shock  dau  to  predict  the  radiation  ema¬ 
nating  from  the  bow  shock. 

Figure  1  depicts  the  major  elements  of  the 
shock  tube  and  instrumenution  used  in  the  experi¬ 
ments,  and  which  are  described  in  deuil  in  Refer¬ 
ence  1. 

Basically,  the  radiation  from  the  incident- 
shock-heated  gas  is  viewed  through  four  windows  at 
the  same  sution  in  the  shock  tube.  Radiometers  in 
the  UV  and  IR  provide  temporally  resolved  proflles. 
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while  the  electronic  spectfometer  provides  UV-vis- 
ible  spectra  recorded  over  selected  time  intervals. 
Wave  speed  and  radiometric  dau  are  simuluneously 
recorded  and  synchronized  via  a  digiul  transient  re- 
corder. 

Examples  of  typical  radiation  profiles  obuined 
with  a  radiometer  at  260  nm  are  presented  in  Fig¬ 
ures  2(a)  and  2(b),  for  the  test  conditions  shown. 
For  all  conditions  of  interest,  the  overall  profiles 
were  similar,  namely,  a  strong,  rapidly  rising  over¬ 
shoot  region,  followed  by  a  decay  to  a  steady-state 
equilibrium.  The  features  of  interest  are  the  u'me- 
to-peak,  peak  signal  levels,  the  decay  rate,  and  the 
signal  levels  at  equilibrium. 

As  will  be  shown  below,  this  overshoot  is  domi¬ 
nated  by  200-320  run  emission  from  the  Gamma 
band  system  of  NO.  The  NO  begins  to  form  at  the 
oitset  of  O2  dissociation  behind  the  shock  front. 
Rapid  NO  electronic  exciution  ukes  place,  owing  to 
the  high  translational  temperature  immediately  be¬ 
hind  the  shock  front.  This  temperature  decays  rap¬ 
idly,  as  it  is  primarily  dominated  by  the  O2 
dissociation  rate.  Thus,  the  electronic  excited  state 
population  is  determined  by  a  number  of  competing 
and  coupled  kinetic  processes,  resulting  in  the  meas¬ 


ured  observables  in  Figure  2.  The  chemical  kinetic 
mechanisms  and  rates  involved  constitute  a  major 
objective  of  this  research. 

In  contrast,  tests  in  pure  oxygen  resiilt  in  pro¬ 
files  such  as  that  of  Figure  2(c).  This  lack  of  over¬ 
shoot  is  significant  to  these  kinetic  studies  and  will 
be  addressed  in  deuil  in  a  later  section. 

As  shown  in  Figure  1,  an  optical  multichannel 
analyzer  (OMA)  was  also  deployed  at  a  window  at 
the  same  shock  tube  stations.  This  iiutrument  could 
be  gated  to  obuin  spectra  over  selected  time  inter¬ 
vals  behind  the  incident  shock  wave.  Typical  spectra 
are  shown  in  Figure  3  where  the  results  of  a  pure 
oxygen  test  are  superimposed  on  that  of  air.  The 
band  structure  of  the  NO  Gamma  system  is  clearly 
identified,  and  can  be  seen  to  override  a  base  of  Oj 
Schumann-Runge  bands. 

By  conducting  a  series  of  tests  in  air  under 
various  initial  shock-tube  pressures  and  shock  ve¬ 
locities,  the  in-band,  volumetric  radiance  of  the 
peak  radiation  was  quantitatively  measured  and  the 
results  reported  in  Reference  1.  The  peak  intensity 
was  found  to  scale  directly  with  altitude  and  was 
shown  to  have  a  very  strong  velocity  dependence. 


Fig.  2  Typical  260  nm  radiance  profilas. 

(a)  air,  Initial  praaaura:  2.25  torr, 
ahock  valocity:  3.81 

(b)  air.  Initial  praaaura:  2.25  torr, 
ahock  valocity:  3.50  mm/^a; 

(c)  oxygon,  initial  praaaura:  2  torr, 
ahock  valocity;  3.89  mm/pa. 


Praaent  Exoarlments 

The  earlier  measurements  described  above 
served  to  establish  the  identification  of  the  species 
and  bands  contributing  to  the  non-equilibrium  ultra¬ 
violet  radiation  from  shock-heated  air,  and  to  meas¬ 
ure  quantitatively  the  dependence  of  this  radiation 
on  altitude  and  shock  velocity.^  The  goals  of  the  re¬ 
cent  studies  primarily  address  the  kinetic  mecha¬ 
nisms  and  rates  that  populate  the  relevant  electronic 
states  and  subsequently  govern  the  radiation  profile. 
These  studies  involved  (a)  the  implementation  of  an¬ 
other  ultraviolet  radiometer  with  higher  spatial  reso¬ 
lution  and  with  enhanced  sensitivity,  (b)  a  series  of 


Fig.  3  Superimpoted  tpoctra  of 

•hock-heatod  air  (upper  trace)  and 
pure  oxygen.  8aa  Raf.  1  for  details. 

measurements  in  which  the  post-shock  densities  of 
NO  and  O2  (and  the  concomitant  temperature  gradi¬ 
ents)  were  varied  over  a  broad  range  for  an  ex¬ 
tended  kinetic  dau  set,  (c)  a  final  analysis  of  the 
equilibrium  NO  and  O2  ra^ation  characteristics,  and 
(d)  the  development  of  a  compuutional  electronic 
sute  excitation  model,  which,  in  conjunction  with 
the  kinetic  dau  set,  aims  to  provide  a  predictive  ca¬ 
pability  for  the  non-equiUbrium  UV  radiation. 
These  program  elements  will  be  addressed  in  the 
balance  of  the  paper. 

The  230  nm  Radiometer 

As  mentioned  above,  the  design  of  this  instru¬ 
ment  sought  to  achieve  higher  spatial  resolution  of 
the  field-of-view  (FOV)  across  the  shock  tube,  as 
well  as  enhanced  sensitivity.  The  design  components 
are  illustrated  in  Figure  4.  The  uniform  FOV  across 
the  tube  is  achieved  by  having  the  aperture  stop  and 
the  image  of  the  field  stop  equal,  at  1  x  8  mm.  At 
this  1  mm  beam  width,  the  O.S  microsecond  digiul 
sampling  time  sets  the  overaU  ideal  resolution  of 
2mm  (at  4  mm/ps).  Recent  analyses  have  indicated 
that  the  overall  resolution  may  be  as  low  as  3-4  mm, 
which  would  be  inadequate  for  high  velocity,  high- 
gradient  shock  radiation  profiles  and  it  it  planned  to 
esublish  this  FOV  in  separate  measurements. 

Another  effect  is  noted,  whereby  Ught  from  the 
advancing  shock,  when  close  to,  but  still  outside,  the 
FOV,  can  scatter  (and  produce  signal)  from  any 
scratch  or  imperfection  in  the  quaru  window  or  lens, 
or  from  the  razor-blade  slit  edges.  This  effect  is 
manifested  in  the  small  pre-shock  "foot*'  teen  in 
most  records. 

Enhanced  sensitivity  was  atuined  by  having  a 
larger  optical  throughput  (cm^-sr)  than  before,  and 
in  using  a  wider  UV  filter  to  capture  more  of  the 
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stronger  NO  Gamma  bands.  Figure  5  presents  the 
relative  filter  functions  of  the  230  nm  and  the  older 
260  run  radiometer.  Comparison  with  the  spectra  of 
Figure  3  iUustrates  the  molecular  bands  subtended 
by  each  of  the  filters. 

Absolute  calibration  was  achieved  in  the  same 
manner  described  in  Reference  1,  in  which  an  NBS- 
traceable,  sundard  deuterium  irradiance  source  is 


deployed.  System  linearity  is  confirmed  by  inverse- 
square-disunce  measurements.  The  maximum  pos¬ 
sible  error  in  the  calibration  is  15%,  although 
self-consistency  tests  with  the  dau  and  ancillary 
measurements  indicate  that  a  probable  error  near 
10%  is  more  appropriate. 

The  Kinetics  Data  Set 

Directly  behind  the  incident  shock  wave  the 
translational  temperature  can  anain  values  near 
6000  K  which  immediately  fall  as  the  Oz  dissociates 
NO  begins  to  form  as  soon  as  O  atoms  become  avail¬ 
able,  and  the  excited  NO  (A)  sute  becomes  popu¬ 
lated  as  it  seeks  to  equilibrate  with  the  decreasing 
temperature.  To  aid  in  elucidating  these  ccnipled 
phenomena,  a  set  of  shock  tube  conditions  was  se¬ 
lected  which  sought  to  separate  and  emphasise  these 
effects.  The  shock  tube  is  ideally  suited  to  such  stud¬ 
ies,  inasmuch  as  the  capability  for  processing  pre¬ 
cisely  controlled  mixtures  can  be  exploited.  In  this 
case,  a  number  of  N2-O2  mixtures  were  used,  with 
initial  pressures  of  S.O.  2.25  and  1.0  torr,  and  with 
shock  velocities  between  3  and  4  mm/|is.  The  largest 
excursion  in  the  post-shock  NO/O2  ratio  was  sought, 
to  provide  the  broadest-range  kinetic  dau  set. 

A  set  of  such  radiation  profiles  measured  with 
the  230  nm  radiometer  is  shown  in  Figure  6,  for  the 
test  conditions  cited  thereon.  The  broad  variation  of 
the  radiation  profiles  anests  to  the  parametric 
changes  among  the  various  tests  and  to  the  useful¬ 
ness  of  the  dau  set  for  studies  of  excited  sute  ki¬ 
netic  populations. 

Also  shown  in  the  figures  are  the  profiles  that 
result  when  the  calculated  O2  Schumaim-Runge 
contribution  has  been  subtracted,  as  described  later. 
Predicting  these  resulunt  NO  Gamma  band  profiles 
is  the  challenge  for  any  scheme  that  addresses  the 
kinetics  of  the  NO  (A)  sute. 

Test  Time 

The  test  time  expected  in  the  shock  tube  can 
be  calculated  from  Mirels’  laminar  theory,  which  ap¬ 
plies  to  all  the  results  reported  herein.  For  the  time- 
hisiories  shown  in  Figure  6,  where  the  initial 
shock-tube  pressure  for  all  cases  was  2.25  torr,  the 
calculated  test  time  decreases  from  145  psec  for  a 
shock  velocity  of  3.24  mm/|isec  to  105  psec  at  3.86 
mm/4sec.  The  observed  test  times,  as  determined 
from  the  termination  of  the  radiation  plateau  in  Fig¬ 
ure  6.  are  somewhat  greater  than  half  this  value, 
which  is  consistent  with  usual  shock-tube  experi¬ 
ence.  The  sensitivity  of  the  UV  radiation  is  such  that 


Fig.  S  Fitter  tranemiaeion  functions. 


IN-BANO  VOLUMETRIC  RADIANCE  (WATTS/CM 


SO  100  150 

(c) 


50  100  150 

(«) 


SO  100  ISO 
(i) 


LABORATORY  TIME 

Fig.  6  Mcaourcd  230  nm  radiation  profilaa  for  thraa  O2  ^N2  mixturat  and  various 
ahock  valoeltias.  Initial  shook  tuba  praasura:  2.25  torr. 

Uppar  eurva:  Obtarvad  radianca  (NO  >♦02  S-R) 

Lowar  eurva;  NOT  contribution  obtalnad  by  subtracting 
ealeulatad  Og  S-R  contribution 


a  drop  in  temperature  of  100  K  decreases  the  radia* 
tion  by  a  factor  of  2. 

4.  Analysis  of  the  Results 


Equilibrium  Ultraviolet  Radiation 

To  calculate  the  contribution  at  wave  number 
w  of  the  equilibrium  spectral  radiation  from  a  given 
(v'.v”)  band,  a  simple  smoothed  Q-branch  structure 
was  used’ 


Iv'v”  (u)  =  3.20  X  10'^°  lR«|^qv'v*  "u 


-Ev'/kT 


(Qv). 


«"ce“Cl“-“<>l 


watts 


cm  ^  sr  cm 
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where 


qv'  V" 

"u 

Ev' 

(Qv), 

c 

Uo 


squared  transition  moment  (atomic 
units) 

Franck-Condon  factor 


number  density  in  the  excited  elec 
tfonic  sute  (molecules/cm’) 


energy  of  vibrational  state  v' 

vibrational  partition  function  of  upper 
electronic  state 

Be'  ) 

Be'  -Be"  *  kT 
bandhead  wavenumber  (cm'l ) 


and  B'  and  B"  are  the  routional  spectroscopic  con¬ 
stants  for  the  upper  and  lower  electronic  sutes.  The 
dependence  of  B  on  v  was  not  included  in  the  calcu¬ 
lation. 

The  total  spectral  radiance  I(u)  was  obtained 
by  summing  these  expressions  over  all  v'.v".  yielding 
spectra  for  the  NO  gamma  and  O2  Schumann-Runge 
systems  as  shown  in  Reference  1.  The  radiation  that 
should  be  measured  by  the  calibrated  detector  was 
calculated  as 

OD 
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where  F(u)  is  the  hlter  transmission  function  and 
D(u)  is  the  detector  response  function. 

Among  the  many  cxperimenul  radiance  pro¬ 
files  obtained  from  shocks  of  varying  velocity  into 
pure  oxygen  and  three  mixtures  of  O2  and  N2.  there 
were  a  number  that  attained  equilibrium  within  the 
available  test  time.  For  each  such  case  Figure  7 
shows  the  equilibrium  value  of  the  cxperimenul  ra¬ 
diance  integral  S  for  the  260  nm  filter,  plotted 
against  the  calculated  value.  Squared  transition  mo- 


CALCULATED  EQUILIBRIUM  RADIANCE 
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Fig.  7  Mtaaurod  and  caleulatad  valuat  of  tha 
260  nm  aquilibrium  radianea  from 
thraa  02/N2  mixturas. 

menu  |ReP  were  chosen  for  the  NO  Gamma  and  O2 
S-R  to  supply  the  best  agreement,  and  are  shown  in 
the  figures.  These  values  are  about  35%  higher  than 
those  reported  in  preliminary  resulu\  and  are  con¬ 
sistent  with  accepted  values  for  the  two  systems.  The 
data  include  resulu  for  the  four  different  02fN2  mix¬ 
tures.  and  a  range  of  shock  speeds  from  3.21  mm/ 
psec  to  4.30  mm/tisec. 

Figure  8  is  a  similar  plot  for  the  230  nm  filter 
data  (available  only  for  the  more  recent  eiqseri- 
menu).  For  this  filter,  the  use  of  a  slightly  different 
set  of  |ReF  values  (as  shown  in  the  figure)  appeared 
to  improve  the  fit. 

The  estimated  precision  in  the  shock-velocity 
measuremenu  is  ±0.2%,  which  corresponds  to  a 
umperature  accuracy  of  about  ±0.2%.  This  can  re¬ 
sult  in  a  variation  in  the  calculated  radiance  of  about 
tS%.  The  ability  to  read  the  ciqierimenul  values  of 
the  radiance  re^u  in  an  unceruinty  of  about  ±3%. 
The  probable  error  in  the  calibration  consunt  is 
±10%. 

Non-Eauilibrium  Ultraviolet  Radiation 

A  time-dependent  record  of  the  230  run  and 
260  nm  radiation  behind  the  shock  wave  was  ob- 
uined  for  four  NilOz  mixtures  (%  O2  *  100.  40. 
22.3.  5)  and  for  a  range  of  velocities  from  3.24  to 
3.86  mm/Rsec.  The  temporal  resolution  of  the  re¬ 
cords  is  estimated  to  be  between  0.5  to  1.0  Rsec. 
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Fig.  8  Maasurad  and  ealeulatad  valuas  of  tha 
230  nm  aquilibrium  radianca  from 
thraa  O2/N2  mixturas. 

To  analyze  the  nonequilibrium  radiation,  it  is 
necessary  to  separate  the  amounts  arising  from  the 
O2  S-R  and  the  NO  Gamma.  The  time-dependent 
experimental  values  of  S  for  the  pure  oxygen  experi¬ 
ments  do  not  show  a  radiation  overshoot  (Figure  2), 
but  rise  quickly  to  steady-state  value.  This  feature 
was  observed  in  work  at  Avco^  in  the  late  1950'$, 
and  attributed  to  the  rapid  predissociation  of  the 
oxygen  excited  B  state.  Ilius  the  B  state  remains  in 
local  thermodynamic  equilibrium  with  the  oxygen 
atom  population,  so  that  the  number  of  O2  mole- 
cules/cm^  in  the  B  state  is  given  by 


where  K  is  the  equilibrium  consunt  for  the  reaction 
and  the  last  factor  in  the  equation  is  an  approximate 
represenution  (appropriate  over  the  range  3000  to 
7000  K)  of  the  required  correction*  to  the  rigid-ro¬ 
tator  simple-harmonic-oscillator  model.  The  O2  ra¬ 
diation  is  taken  to  be  proportional  to  RgjfB) ' 
determined  from  the  equilibrium  measurements. 


We  have  assumed  that  the  O2  S-R  system  fol¬ 
lows  the  same  behavior  in  the  radiation  from  N2fQ2 
mixtures,  with  the  Q2(B)  sute  suying  in  local  equi¬ 
librium  with  the  oxygen  atom  concentration  at  the 
translational  temperature.  The  oxygen  atom  concen¬ 
tration  behind  the  shock  was  calculated  from  a 
normal-shock  nonequilibrium  code*'\  and  this 
population  was  used  to  determine  the  B  state  partici¬ 
pation  in  the  observed  radiation.  This  amount  of 
radiation  was  subuacted  from  the  toul  observed, 
and  the  remainder  was  attributed  to  the  NO  Gamma 
system. 

In  the  results  shown  in  Figure  6,  the  total  ra¬ 
diation  and  the  portion  attributed  to  NO  Gamma  are 
shown.  These  graphs  of  the  NO  radiation  provide  a 
record  of  the  population  of  the  NO  (A)  sute  as  a 
function  of  time. 

The  NO  (A)  sute  could  be  produced  either  di¬ 
rectly  or  by  transfer  from  the  meusuble  NafA) 
sute.  Several  reporu*'^  have  speculated  on  the  pos¬ 
sible  candidate  reactions  for  both  NO  (A)  and  N2(A) 
production.  The  reactions  that  are  considered  likely 
for  direct  production  of  the  NO  (A)  sute  are  the 
electron  collisional  exciution  of  NO(X)  and  the  re¬ 
combination  of  N  and  O  into  the  NO(A)  sute.  The 
elecuon  density  in  the  present  eiqieriments  is  very 
low,  so  that  the  first  reaction  cannot  be  considered 
as  likely.  The  recombination  reaction  does  not  pro¬ 
duce  an  overshoot,  and  so  will  not  supply  the  ob¬ 
served  profiles.  Heavy-body  collisional  exciution, 
either  of  NO(A)  directly  or  of  N2(A)  followed  by 
transfer  to  NO  (A),  will  give  overshoots  of  the  type 
that  have  been  observed  in  these  experiments.  How¬ 
ever,  no  definitive  set  of  rates  that  will  describe  the 
experimenul  results  shown  in  Figure  6  has  yet  been 
determined. 

We  have  begun  the  design  of  an  LIF  technique 
to  monitor  the  production  of  the  meusuble  N2(A) 
sute.  Near-IR  absorption  into  the  N2(B)  sute  with 
observation  of  subsequent  First  Positive  band  emis¬ 
sion  is  planned.  A  direct  comparison  of  the  N2(A) 
and  NO(A)  sute  profiles  will  provide  a  significant, 
complemenury  aid  in  studying  these  complex  exciu¬ 
tion  kinetics. 

Equilibrium  Infrared  Radiation 

The  infrared  radiation  from  the  shock-heated 
gas  was  measured  using  a  HgCdTe  detector  and  a 
narrow -bandpass  filter  with  consunt  transmission 
from  1749  cm**  to  1954  cm**,  allowing  good  re¬ 
sponse  for  the  NO  fundamenul  IR  band,  which  is 
centered  at  1850  cm**.  The  sensitivity  of  the  detec¬ 
tor  was  assumed  consunt  over  this  spectral  range. 
The  measured  signals  contained  an  apipreciable 
amount  of  noise,  but  for  many  conditions  a  plateau 


region  was  reached  where  a  final  signal  value  could 
be  determined.  An  example  is  shown  in  Figure  9. 
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Fig.  9  Measured  infrared  rediation  with 
equiiibrium  piateau. 

Initial  shock  tube  pressure;  5.00  torr 
Shock  velocity:  3.42  mm/^s. 

The  strength  of  the  signal  to  be  expected  can 
be  obtained  by  calculating  the  radiance  of  a  single. 
v'.J';  v",J"  transition,  and  summing  these  values 
over  all  the  lines  within  the  filter  bandpass.  The  re¬ 
sulting  signal  suength  is  less  dependent  on  tempera¬ 
ture  than  is  the  toul  band  radiance,  because  at 
higher  temperatures  a  larger  portion  of  the  band  is 
outside  the  range  of  the  filter.  This  calculated  signal 
suength  can  be  written  as: 

Sir  =  ^  (T-3000) 

^  -  2.5  X  10’* (T-3000)^}  (4) 

over  the  temperature  range  3000  to  5000  K,  where 

X  B  pathlength  s  7.62  cm 
n^  *  2.689  X  10  ’’ molecules/cm  ^ 

and  is  the  number  density  of  NO  molecules. 

A  is  a  calibration  constant  for  the  system  which  was 
not  determined.  However,  the  consistency  of  the 
data  can  be  shown  by  matching  the  calculated  results 
to  a  single  datum  point  and  using  this  coiuunt  for 
the  other  points.  The  resulting  plot  for  all  erq^eri- 
ments  that  clearly  reached  an  equilibrium  plateau  is 
shown  in  Figure  10.  The  variations  shown  are  consis¬ 
tent  with  the  noise  level  of  the  infrared  signal. 
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Fig.  10  Maaturad  and  calculatad  valuaa  of  tha 
aquillbrium  Infrarad  radiation  from  two 
O2/N2  mixturas. 

Non-Etfuilihrium  Infrared  Radiation 

Three  examples  of  time-dependent  infrared 
signals  are  shown  in  Figure  1 1 .  These  examples  are 
from  the  runs  for  which  the  230  nm  radiation  is 
shown  in  Figure  6.  The  expected  signal  is  shown  by 
the  smooth  curves,  which  are  calculated  using  equa¬ 
tion  (4)  and  values  of  ni,jQ  and  T  determined  by 
the  normal-shock  calculation^  ',  scaled  to  the  equi- 
Ubrium  value  (Figure  10). 

The  lack  of  spatial  resolution  in  the  infrared 
measurements  is  demonstrated  by  the  inability  of  the 
system  to  follow  the  calculated  radiation  overshoot, 
and  by  the  sizable  “foot*'  generated  before  arrival  of 
the  shock  wave  at  t  >  0.  With  this  poor  resolution  it 
would  also  be  expected  that  internal  reflections  in 
the  shock  tube  might  contribute  to  the  resolution 
problem.  However,  further  narrowing  of  the  optical 
beam  led  to  an  unacceptable  increase  in  signal  noise. 

To  address  this  problem,  a  new  infrared  radi¬ 
ometer  has  been  designed  and  insulled  in  the  shock 
tube.  It  utilizes  an  indium  antimonide  detector, 
which,  despite  the  response  roll-off  at  5.5  pm.  still 
provides  a  factor  of  10  gain  in  overall  SNR  over  that 
of  the  HgCdTe  detector  deployed  previously.  In  ad¬ 
dition,  the  spatial  resolution  has  been  enhanced  by 
means  of  appropriate  slits  and  CaF2  lenses.  An  ex¬ 
tended  range  of  IR  measurements  at  higher  SIVR  is 
anticipated  for  future  tests. 
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^—Summary 


Calculated  predictions  of  non-equilibrium  ul¬ 
traviolet  radiation  from  the  bow  shock  of  hypersonic 
vehicles  depends  critically  on  relevant  e]q>erimenul 
data,  srfiich  comprise  the  major  goal  of  these  studies. 
Earlier  work  led  to  the  identifrcation  of  the  radiating 
species  and  to  the  velocity  and  altitude  dependence 
of  the  UV  radiances. 

Further  experimenution  is  reported  that  ad¬ 
dresses  the  kinetics  of  the  relevant  electronic  sutes. 
In  particular,  an  extensive  kinetic  data  set  has  been 
obtained,  over  which  the  NO  and  O2  post-shock 
densities  have  been  broadly  varied.  A  compuutional 
model  has  been  developed,  involving  the  reactions 
and  rates  for  excited  sute  populations.  Iterative 
comparative  efforu  are  now  underway  to  reconcile 
the  calculated  and  measured  nonequUibrium  radia¬ 
tion  profiles. 

Equilibrated  radiance  levels  for  some  cases 
permitted  an  evaluation  of  the  NO  Gamma  and  O2 
Schumann-Runge  uansition  moments,  from  which 
their  relative  contributions  to  air  radiation  were  de- 
urmined.  Similarly,  equilibrium  infrared  measure¬ 
ments  were  used  to  esublish  that  the  NO 
fundamental  band  sequence  provided  a  quantiutive 
measure  of  the  NO  concentration.  Additional  and 
improved  diagnostics  are  planned  to  record  further 
observables  for  the  study  of  these  kinetic,  radiative 
processes. 


Fig.  11  Maaaurad  and  ealeulatad 

non-aquillbrium  Infrarad  radiation 
profilas  for  thraa  ahoek  valoeltias  in 
40%  02/60%  N2  mixturas. 

Initial  ahoek  tuba  prastura:  2.25  torr 
Taat  tima  ands  at  -70  us.  (Saa  Fig.  6) 
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APPENDIX  D 
IR  SPECTRAL  ANALYSIS 


This  appendix  presents  a  summary  of  the  equations  that  were  used  in  the  analysis  on  the 
equilibrium  IR  radiation  from  nitric  oxide,  and  applies  them  to  the  results  obtained  in  the  current 
experimc  The  equations  are  discussed  in  detail  in  References  7  and  13. 


The  wave  numbers  of  the  rotational-vibrational  lines  in  the  P,  Q,  and  R  branches  are  given 


D-1)  =1^,  -2.('tV()B^ -oC^[(rj+(f-2(jV0v''  ~ij 

D-2)  Vq  = 

♦ 

D-3)  2.jV  ~\) 

for  a  transition  from  the  state.  The  wavenumber  l)^i  is  given  by 


J=l,2,3,  .  .  . 


J=l,2,3,  .  .  . 


J=2,3,4,  . 


D-4)  •=  CJe-2-^^COe'Xg  -i-3.25>  =[jqo^i.(P3  -  v'-(o.o03>9)  cm 

The  constants  for  the  NO  molecule  are; 

Be  =  1.7046  cm"'  and  ^TCe  =  0.0178  cm'* 

The  intensities  of  the  lines  are  given  by  - 

,  _  (jv- 

T*  _  f Cer>i  V  ^  k,T 

D-5)  -e- 

n./  yi  ■.  Q  j/'  -  Cy'+i)!  __  _<3v_' 

n-«  -e 

D-7)  -f-R -  qTc;  ®  ^ 

Watt/ H0L€:cui.e-£^ 

WHEKE  (Sv=^[«e-(''^f')<:.7e'’fe*(''^+i.5v+o.7s-;aJeUeJ 


and  k=0.fc95  cnr^K  ^  Cem  is  the  emission  constant  for  the  NO  band  system  and  a  value  of 
Cem“2.97±0.54  x  watts  cmVmolecule  -5»*  was  reported  in  References  7  and  13. 


The  partition  functions  are  approxmiated  as 
D-8)  Qr  =  - — — 7*7——. 

and  ,  V. 

D-9)  Qv  =  (  ( -  e”  feT  ) 

The  total  in-band  radiation  measured  is  then  calculated  as 

D-IO)  R  =  '>^MolZ[lpfVf>*-Ie'v«*lR'VR] 

where  Fv  is  the  combined  filter-detector  response  at  wavenumber  i)  ,  as  described  in  Section  3.2, 
■n-ivjo  is  the  number  density  of  NO  molecules  (molecules/cm^)  and  Jis  the  diameter  of  the  shock 
tube,  7.62  cm.  The  values  of  Fv  were  taken  as  constants  over  specified  wavenumber  intervals,  as 
given  in  Table  D- 1. 


Wavenumber  Interval  (cm*i) 

Average  Response 
Fy 

1724.1 

1773.0 

0.00445 

1773.0 

1793.7 

0.0311 

1793.7 

1816.5 

0.11565 

1816.5 

1839.9 

0.2975 

1839.9 

1863.9 

0.515 

1863.9 

1888.6 

0.7105 

1888.6 

1913.9 

0.8175 

1913.9 

1938.0 

0.851 

1938.0 

1964.6 

0.8815 

1964.6 

1992.0 

0.7615 

1992.0 

2020.2 

0.395 

2020.2 

2049.2 

0.08fc4 

2049.2 

2079.0 

0.0089 

TABLE  D- 1 .  FILTER-DETECTOR  RESPONSE 


The  values  of  Fv  are  relative  numbers,  and  are  the  same  values  used  in  the  calibration  of  the 
system.  As  shown  in  Figure  4,  this  calibration  yields  the  relation 

-4- 

D-11)  S-  1.04-><10  V"  W/oi*--5R 

where  S  is  the  in-band  radiation,  /  IbbF^^.  and  V  is  the  number  of  millivolts  measured  by  the 
detector. 

The  values  obtained  for  the  calculation  indicated  in  Eqn.  D-10  are  shown  in  Table  D-2, 
calculated  for  three  temperatures  covering  the  range  of  interest  in  the  present  experiments. 


T(“K) 

Cem 

3000 

3.507  X  10*2 

4000 

4.755  X  10*2 

5000 

5.758  X  10*2 

TABLE  D-2.  CALCULATED  RADIATION  INTEGRALS 
These  values  of  R/n^CCem  can  be  expressed  as 


D-12)  _ i _  =  'f(T)=  3*5'07>^l0  Fl-i-SiZIxlO  (T-BDeo)! 

Turning  now  to  the  experimental  results  for  the  infrared  measurements  with  NO,  the  value 
of  S  for  each  experiment  is  obtained  from  the  number  of  millivolts  measured  at  equilibrium,  by 
using  Equation  D-11.  Using  i,  *  7.62  cm,  the  value  of  S/Jt  [1+3.21  x  lO-^T-3000)]  is  plotted  in 
Figure  D- 1  vs.  dno*  According  to  Equation  D-12  this  should  yield  a  straight  line  with  slope  equal 
to  3.506  X  10*2  Qjjj  From  Figure  D- 1  a  value  of  Cem  is  obtained  as 

Cem  “  3.66  X  10-33  watts  cm'^/molecule  -sr 
slightly  higher  than  the  value  quoted  in  References  7  and  13. 
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Full-Scale  Records  of  Data  in  Figure  14 


-20 


X10»**-9  X10» 


TORY  Tire*  rucROSEca DS 


X10***»8  X10**^*»8 


